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Abstract. This paper concerns the long-term behavior of population systems, and in particu- 
lar of chemical reaction systems, modeled by deterministic mass-action kinetics. We approach two 
important open problems in the field of Chemical Reaction Network Theory, the Persistence Con- 
jecture and the Global Attractor Conjecture. We study the persistence of a large class of networks 
called lower- endotactic and in particular, we show that in weakly reversible mass-action systems with 
two-dimensional stoichiometric subspace all bounded trajectories are persistent. Moreover, we use 
these ideas to show that the Global Attractor Conjecture is true for systems with three-dimensional 
stoichiometric subspace. 
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1. Introduction. Mass-action systems are a large class of nonlinear differential 
equations, widely used in the modeling of interaction networks in chemistry, biology 
and engineering. Due to the high complexity of dynamical systems arising from non- 
linear interactions, it is very difficult, if not impossible, to create general mathematical 
criteria about qualitative properties of such systems, like existence of positive equi- 
libria, stability properties of equilibria or persistence (non-extinction) of variables. 
However, a fertile theory that answers this type of questions for mass-action systems 
has been developed over the last 40 years in the context of chemical reaction systems. 
Generally termed Chemical Reaction Network Theory [HJII1II1II1II3IISII1I1I] 
this field of research originated with the seminal work of Fritz Horn, Roy Jackson and 
Martin Feinberg [III US 120] and describes the surprisingly stable dynamic behavior 
of large classes of mass-action systems, independently of the values of the parameters 
present in the system. This fact is very relevant, since the exact values of the system 
parameters are typically unknown in practical applications. Although the results in 
this paper will be applicable to general population systems driven by mass-action ki- 
netics, they will be developed within the frame of Chemical Reaction Network Theory. 

A large part of this paper is devoted to persistence properties of mass-action 
systems. A dynamical system on R™ is called persistent if forward trajectories that 
start in the interior of the positive orthant do not approach the boundary of R> (see 
section 2.6 for a rigorous definition). Note that, throughout this paper, trajectory 



will always mean bounded trajectory. For systems with bounded trajectories, this is 
equivalent to saying that no trajectories with positive initial condition have w-limit 
points on the boundary of R> . Persistence answers important questions regarding 
dynamic properties of biochemical systems, ecosystems, or infectious diseases, e.g. 
will each chemical species be available at all future times; or will a species become 
extinct in an ecosystem; or will an infection die off. One of the major open questions 
of Chemical Reaction Network Theory is the following: 
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Persistence Conjecture [10 . Any weakly reversible mass-action system is persis- 
tent. 

A weakly reversible mass-action system is one for which its directed reaction graph 
has strongly connected components (Definition 2.3 1. A version of this conjecture was 
first mentioned by Feinberg in [131 Remark 6.1.E]; that version only requires that no 
trajectory with positive initial condition converges to a boundary point. A stronger 
version of the Persistence Conjecture (called the Extended Persistence Conjecture) was 
formulated by Craciun, Nazarov and Pantea in |10j and it was shown to be true for 
two-species systems. Moreover, in that case, weakly reversible mass-action systems 
are not only persistent, but also permanent (all trajectories originating in the interior 
of R™ eventually enter a fixed compact subset of the interior of K> )- 

In recent approaches to the Persistence Conjecture, the behavior of weakly re- 
versible mass-action systems near the faces of their stoichiometric compatibility classes 
(minimal linear invariant subsets) was considered. It is known that w-limit points may 
only lie on faces of the stoichiometric compatibility class that are associated with a 
semilocking set pQ (see also [131 Remark 6.1.E]), or siphon in the Petri net litera- 
ture [5J [57] . Anderson [I] and Craciun, Dickenstein, Shiu and Sturmfels [9J showed 
that vertices of the stoichiometric compatibility class cannot be w-limit points. More- 
over, Anderson and Shiu [4] proved that for a weakly reversible mass-action system, 
the trajectories are, in some sense, repelled away from codimension-one faces of the 
stoichiometric compatibility class. 

In this paper we prove the following version of the Persistence Conjecture for 
systems with two-dimensional stoichiometric compatibility classes (Theorem 5.1): 
Theorem 5.1. Any n-variable mass-action system with bounded trajectories, two- 
dimensional stoichiometric compatibility classes and lower- endotactic stoichiometric 
subnetworks is persistent. 

Here a stoichiometric subnetwork is a union of connected components of the reac- 



tion graph (see Definition 2.8) and the requirement of lower- endotactic (Definition 



3.4) stoichiometric subnetworks is less restrictive than that of weak reversibility. We 



suggest that the hypothesis of "lower-endotactic" arises naturally in the context of 
persistence of mass-action systems. Moreover, n-variable mass-action is a general- 
ization of mass-action where each reaction rate parameter is allowed to vary within 
a compact subset of (0, oo) (see Definition 2.6). Therefore this theorem implies the 
following: 

Corollary. Any weakly reversible mass-action system with two-dimensional stoichio- 
metric compatibility classes and bounded trajectories is persistent. 
Note that our proof of Theorem 5.1. above (and of its corollary) requires the hypoth- 
esis of bounded trajectories. However, it has been conjectured that all trajectories 
of weakly reversible mass-action systems are bounded [3], and this conjecture has 
recently been proved for networks whose reaction graph has a single connected com- 
ponent [3]. Also, a stronger statement is known to be true for two-species networks: 
any endotactic, n-variable mass-action system with two species has bounded trajec- 
tories [TP] . 

The Persistence Conjecture is strongly related to another conjecture which is of- 
ten considered the most important open problem in the field of Chemical Reaction 
Network Theory [2 El |9l [10], namely the Global Attractor Conjecture. This con- 
jecture was first formulated by Horn [19] and concerns the long-term behavior of 
complex-balanced systems, i.e., systems that admit a positive complex-balanced equi- 
librium (see Definition 2.10). Horn and Jackson showed that if a mass-action system 
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is complex-balanced then there exists a unique positive equilibrium in each stoichio- 
metric compatibility class and this equilibrium is complex-balanced |20j . Moreover, 
each such equilibrium is locally asymptotically stable in its stoichiometric compati- 
bility class due to the existence of a strict Lyapunov function |20j . In two subsequent 
papers [Til EH] Feinberg and Horn showed that weakly reversible mass-action sys- 
tems which are also deficiency zero (Definition 2.11) are complex-balanced. This 
fact is remarkable since it reveals a wide class of reaction systems that are complex- 
balanced only because of their structure and regardless of parameter values. For a 
self-contained treatment of Chemical Reaction Network Theory, including the results 
mentioned above, the reader is referred to |12j . 

The Lyapunov function of Horn and Jackson does not guarantee global stability 
for a positive equilibrium relative to the interior of its compatibility class. This fact 
is the object of the Global Attractor Conjecture: 

Global Attractor Conjecture. In a complex-balanced mass-action system, the 
unique positive equilibrium of a stoichiometric compatibility class is a global attractor 
of the interior of that class. 

It is known |12j that complex-balanced systems are necessarily weakly reversible. 
On the other hand, trajectories of complex-balanced systems converge to the set of 
equilibria [I, 24, 26J, so it follows that the Persistence Conjecture implies the Global 
Attractor Conjecture. 

A series of partial results towards a proof of the Global Attractor Conjecture 
have been obtained in recent years. It is known that the conjecture is true for sys- 
tems with two-dimensional stoichiometric compatibility classes ([4]; see also the recent 
work of Siegel and Johnston [IS]), and for three-species systems (TUj. Recently, An- 
derson proved that the conjecture holds if the reaction graph has a single connected 
component [2]. 

In this paper we prove the Global Attractor Conjecture for systems with three- 
dimensional stoichiometric compatibility classes (Theorem 6.3): 

Theorem 6.3. Consider a complex-balanced weakly reversible mass-action system 
having stoichiometric compatibility classes of dimension three. Then, for any positive 
initial condition Co, the solution c{t) converges to the unique positive equilibrium which 
is stoichiometrically compatible with cq. 

Aside from being significant in the field of polynomial dynamical systems and 
relevant in important biological models [T3J [Tj)J [Ml US] , Chemical Reaction Network 
Theory, and in particular the two conjectures discussed above, have ramifications in 
other well-established areas of mathematics. For example, [9] stresses the connec- 
tion with toric geometry and computational algebra; in that work complex-balanced 
systems are called toric dynamical systems to emphasize their intrinsic algebraic struc- 
ture. Also, [33J studies the rich algebraic structure of biochemical reaction systems 
with toric steady states. Furthermore, the unique positive equilibrium in a stoichio- 
metric compatibility class of a complex-balanced system is sometimes called the Birch 
point in relation to Birch's Theorem from algebraic statistics [9l [2Tj. 

This paper is organized as follows. After a preliminary section of terminology 
and notation, we introduce the lower-endotactic networks in section [3] and follow 
with a discussion of our main technical tool, the 2D-reduced mass-action system in 
section [4j The main persistence result of the paper is contained in section [5] (Theorem 



5.1) and our result on the Global Attractor Conjecture (Theorem |6.3[) is proved in 



section [6] A critical part in the proof of the latter theorem resides in the result of 

3 



Theorem |6.2| which analyzes the behavior of weakly reversible mass-action systems 
near codimension-two faces of stoichiometric compatibility classes. 



2. Preliminaries. A chemical reaction network is usually given by a finite list 
of reactions that involve a finite set of chemical species. An example with four species 



A,B,C and D and five reactions is given in (2.1 1 



B + D^A + C^A + B^C + D 2A^A + D^2D (2.1) 

The interpretation of the reaction A + B ^> C + D, for instance, is that one molecule 
of species A combines with one molecule of species B to produce one molecule of 
each of the species C and D. The objects on both sides of a reaction are formal 
linear combinations of species and are called complexes. According to the direction 
of the reaction arrow, a complex is either source or target. This way, the reaction 



A + B — » C + D in (2.1) has A + B as source complex and C + D as target complex. 
The concentrations ca,cb,cc and cd vary with time by means of a set of ordinary 
differential equations, which we will explain shortly. In this preliminary section of 
the paper we review the standard concepts of Chemical Reaction Network Theory 
(see [12] ) and introduce some new terminology that will be useful further on. In 
what follows, the set of nonnegative, respectively strictly positive real numbers are 
denoted by ]R>o and ]R>o . For an integer n > 1 we call ]R™ the positive orthant. 
The boundary of a set K C W 1 will be denoted by dK and the convex hull of K will 
be denoted by conv(K). Also, we will denote the transpose of a matrix A by A*. 

2.1. Reaction networks. If 7 is a finite set then we denote by Z> and R> 

the set of all formal sums a = a,i where oti are nonnegative integers, respectively 

iei 

nonnegative reals. 

Definition 2.1. A chemical reaction network is a triple (S,C,1Z), where S is 
the set of species, C C Z> is the set of complexes, and 1Z is a relation on C, denoted 
P — » P' , representing the set of reactions of the network. The reaction set 1Z cannot 
contain elements of the form P — > P and each complex in C is required to appear in 
at least one reaction. 

For simplicity, we will often denote a reaction network by a single letter, for 
instance M = (S,C,1Z). For technical reasons we have chosen to neglect a third 
requirement that is usually included in the definition of a reaction network (see [12 ): 
each species appears in at least one complex. This condition is not essential in the 
setting of this paper. 



In (2.1) the set of species is S = {A, B,C, D}, and the set of complexes is C = 
{B + D,A + C,A + B,C + D,2A,A + D, 2D}. 

Once we fix an order among the species, any complex may be viewed as a column 
vector of dimension equal to the number of elements of S. For example, the complexes 
A + B and 2D in (O) may be represented by the vectors (110 0)*, and (0 2)*. 



With this identification in place, we may now define the reaction vector of a reaction 
P -> P' e TZ to be P' - P. 

Definition 2.2. The stoichiometric subspace of the reaction network JV = 
{S,C,K) isS= span{P' - P | P -> P' e 11}. 



Example 1 . The stoichiometric subspace of the reaction network (2.1) is the 
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column space of the stoichiometric matrix 



A = 



( 1 -1 0-1-1 1 \ 
-111-10 
1-1-110 

\ -1 1 1 1 -1 J 



It is easy to see that S = {(a, b, —b, — a)*|(a, b) £ M 2 }. 

Any reaction network Af can be viewed as a directed graph whose vertices are 
the complexes of Af and whose edges correspond to reactions of Af. Each connected 
component of this graph is called a linkage class of Af. 

Definition 2.3. A reaction network Af is called weakly reversible if its associated 
directed graph has strongly connected components. 

In other words, Af is weakly reversible if whenever there exists a directed arrow 
pathway (consisting of one or more reaction arrows) from one complex to another, 
there also exists a directed arrow pathway from the second complex back to the first. 

2.2. Reaction systems. Throughout this paper we let n denote the number of 
species of a reaction network Af — (S,C,TZ), we fix an order among the species, and 
we denote S — {A 1; . . . , X n }. We also let c(t) eR 5 = R™ denote the (column) vector 
of species concentrations at time t > 0. From here on, "vector" or "point of K n " will 
always mean "column vector" , even if, for simplicity, the notation * may not always 
be used. The concentration vector c(t) is governed by a set of ordinary differential 
equations that involve a reaction rate function for each reaction in 1Z. 

Definition 2.4. A (non-autonomous) reaction system is a quadruple (S, C, 1Z, K) 
where Af — (S,C,7Z) is a reaction network with n species and K : IR> x M" — > R> 
is a piecewise differentiable function called the kinetics of the system. The component 
Kp^pi of K is called the rate function of reaction P — > P'. Letting P = (mi, . . . , m n ) 
and x = (xx, . . . , x„), Kp^p> is assumed to satisfy, for any t > 0, the following: if 
Xi = and wii ^ then Kp^pi(t,n) = 0. The dynamics of the system is given by the 
following system of differential equations for the concentration vector c(t) : 

c(t)= Kp^ P ,(t,c(t))(P' -P). (2.2) 
p^P' 

Note that we will often use the short notation (Af, K) for a reaction system. 
The regularity condition on K may be replaced by any other condition that 



guarantees uniqueness of solutions for (2.2). Sometimes, additional properties are 



required of K [31 El E] ■ For example, it is commonly assumed that if the ith species is 
not a reactant in P — > P' (i.e. m, = 0) then Kp^p* does not depend on xi. Another 
widespread assumption is that K is increasing with respect to reactant concentrations, 
i.e. -r^r-Kp^pi > if nrii ^ 0. These conditions are automatically satisfied for the 
kinetics treated in this paper. 
If c e M> we let 

T(c ) = {c(t) | t > 0,c(0) = c } 

denote the trajectory of (Af,K) with initial condition cq. If K(t,x.) = K(x) does not 
depend explicitly on time, we say that c* G IR" i s an equilibrium of the reaction 



system (Af,K) if it is an equilibrium of the corresponding differential equations (2.2) 



Note that the condition on Kp^pi imposed in Definition |2.4| makes the nonneg- 



ative orthant M™ forward invariant for (2.2). Under mild additional assumptions on 



K, the positive orthant M" is also forward-invariant for (2.2) (see [26]). For exam- 
ple, this will be the case for n-variable mass-action kinetics, the main type of kinetics 
considered in this paper (Definition 2.6). 
Integrating (2.2) yields 

c(f)=c(0)+ J2 ( [ K P ^ p/ (s,c(s))ds)(P' -P) 
p^P' ' 

and it follows that c(t) is contained in the affine subspace c(0) + S for all t > 0. 
Combining this with the preceding observation we see that (c(0) + S) nK™ is forward 
invariant for (2.2). 

Definition 2.5. Let c e R n . The polyhedron (c + S) D M> is cashed i/ie 
stoichiometric compatibility class of Cq. 

Note that, throughout this paper, "polyhedron" will always mean "convex poly- 
hedron", i.e., an intersection of finitely many half-spaces. To prepare for the next 
definition we introduce the following notation: given two vectors u,v € , we 

n 

denote u v = J^J , with the convention 0° = 1. 

Definition 2.6 ([Hi])- A k- variable mass-action system is a reaction system 
(J\f,K) where J\f — (S,C,1Z) and the rate function of P — > P ElZ is given by 

Kp^ P ,(t,x) = K P ^p,(t)x p . (2.3) 

Here K : R>o — > (t)> l/ 7 ?) 7 ^ f or some rj < 1 is a piecewise differ entiable function called 
the rate-constant function. 

To emphasize the rate-constant function, we denote a K-variable mass-action sys- 
tem by (S,C,TZ, k), or by (Af, k). Note that if the rate-constant function is fixed in 
time, K-variable mass-action becomes the usual mass-action. A few biological ex- 
amples of K-variable mass-action models that are not mass-action are presented in 

Therefore, a n- variable mass-action system gives rise to the following non-autonomous, 
and usually nonlinear, system of coupled differential equations: 

c(t)= Kp^P'(t)c(t) p (P' -P). (2.4) 
p->P' 



Example 2. We endow the reaction network (2.1) with n-variable mass-action 
kinetics of rate- constant function specified on the reaction arrows in {2.5). 

B + D^A + C^A + B^C + D 2A A + D ^ 2D (2.5) 

Ka(t) 

We have c(t) = (c/i(i), cb(£), Cc(t) , cpj(t)) and note that, for example, c(t) A+B = 
c(£) (1 l o o) = CA (t)c B (t). From (24) we have 

C = Ki(t)cBCDA 1 + K2(t)cACcA2 + K3(t)cACcA 3 + K i (t)c A CBA i + K 5 (t)c\A 5 + Ke(t)c 2 D A 6 

for all t > 0, where A4 is column i of the stoichiometric matrix A given in Q), i.e., 
the reaction vector of reaction i. Therefore the differential equations corresponding to 



(2.5) are 



CA = H\{t)c B CD ~ K 2 (t)c A C C - Ki(t)c A CB ~ K 5 (t)c A + K 6 (t)c D (2.6) 

C-B = -K\(t)c B CD + K 2 (t)c A C C + K 3 (t)c A C C ~ H^(t)c A C B 

CC = Kl(t)c B C D - K 2 (t)c A C C ~ K 3 (t)c A C C + K 4 (t)c A C B 

CD = -Kl(t)c B CD + K 2 (t)c A C C + K±(t)c A CB + H5(t)c A - K 6 (t)c D , 



2.3. Sums of reaction systems. A reaction network (S,C,7Z) is called a sub- 
network of (£>,' C ,1Z') if S C S 1 , C C C' and 1Z C TZ 1 . H1Z' has an associated kinetics 
K then restricting K to reactions of TZ defines a kinetics for TZ. On the other hand, 
if 7V" S = (<S S , C s , 7?. s ) , s G {1, ...,p} are reaction networks, their union, denoted by 
U^=i(>5 s ,C s ,7^ s ) or simply by [X =1 Af s , and defined as the triple (Ug=i $s, Us=i C s, Us=i 
is also a reaction network. If each Af s has an associated kinetics K s , we can define a 
kinetics for Us=i-^s by simply adding all K s . 

Definition 2.7. The sum of the reaction systems (S S ,C S ,7Z S , K s ) is the reaction 
system (S,C,TZ,K) where (S,C,TZ) = {J P S=1 {S S ,C S ,TZ S ) and 

Kp^pi (t, x) = K S: p^ P ,(t,x) 

for and all (f,x) eR> xR£ . We will denote this (S ,C,1Z, K) by\J p s=1 (S s ,C s ,K s ,K s ) 
or simply by {J P S=1 {N S , K s ), where N s = (S S ,C S ,TZ S ). 

For example, any reaction system is the sum of the reaction systems corresponding 
to its linkage classes. Similarly, any reaction system is the sum of the reaction systems 
corresponding to its stoichiometric subnetworks, which we define next. 

Definition 2.8. A reaction network (S,C,TZ) with stoichiometric subspace S 
can be written uniquely as a union of subnetworks 

p 

(S,C,1Z) = \J(S S ,C S ,1Z S ) 

s=l 

where {C s } s e{i,...,p} is a partition of C such that two complexes in C are in the same 
block of the partition if and only if their difference is in S. We call each (S S ,C S ,7Z S ) 
a stoichiometric subnetwork of (S,C,1Z). 

Example 3. The diagram in Figure \2~1\ represents the reaction network 

C -> B -> A, 2B ^ A + B, B + C -> 2A -> 2C. 

This reaction network has three linkage classes, and two stoichiometric subnetworks 

{C -> B -> A} and {IB ^ A + B, B + C ->• 2A -+ 2C}. 



Note that there exist vectors a\, . . . ,a p <E R" such that for all s e {1, . . . ,p}, we 
have C s C a s + S and the affinc subspaccs a s + S, s £ {1, . . . ,p} are pairwise disjoint. 
Also note that each stoichiometric subnetwork is a union of linkage classes. 

Example 4. The reaction network (2.1) has two linkage classes which belong to 
two different stoichiometric subnetworks since (A+B) — 2A = (—1, 1, 0, 0) ^ S (recall 
S from Example [IJ). Therefore the reaction network has exactly two stoichiometric 
subnetworks which coincide with its linkage classes. 
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2.4. Projected reaction systems. For W C {1, . . . , n} we define 

n w : E" -> M w 

to be the projection onto W, i.e. the orthogonal projection onto R w . 

Definition 2.9 ([22]; see also 0). LetAf = (S,C,Tl) be a reaction network with 
S = {Xi,...,X n } and let W C {1, . . . ,n}. Let S w = {X, | i e W}, 

Tl w = {tt w (P) ->■ n w (P') \P-> P' EK, such that n w (P) ^ tt w (P')} 

andCw C ttw(C) be the set of complexes inlZw- The reaction network (Sw ,Cw ,T?-w)i 
is called Af projected onto W and denoted n-\y(Af). 

In other words, the projection of (S,C,1Z) onto W is obtained by deleting the 
species X^i € CW from all reactions in 1Z and further removing the resulting reactions 
for which the source and the target complexes are the same. Here and from now on 
ZW denotes the complement of W in {1, ... , n}. 

If (Af, K) = (S,C, 1Z, K) is a reaction system and c(t) = (c\(t), . . . , c n (t)) is a so- 



lution of the corresponding system of differential equations ( 2.2 1 with initial condition 
Co € K." , then ww(c)(t) is a solution of the following system of differential equations: 

j*w(c)= £ ( ]T K P ^ P ,(t,^ w (c))\(Q'-Q) (2.7) 

■n w {P)=Q.-n w {P')=Q'} 



with initial condition Ttwipo)- Equation \2.7\ is obtained from (2.2 1 by (i) keeping 
only the equations for Cj with i € W; (ii) writing Kp^pi(t,c) = K p^p>(t,Trw(c)) to 
illustrate that Cj, i € CW, are written either in terms of Cj, i <E W, or as functions of 
£; and (m) lumping together the rates of reactions P — > P' that project to the same 
reaction in TZw- The system of differential equations (2.7) defines a kinetics K\y 



for iiw(Af), where, for any reaction Q — > Q' g 72.^/, if^Q^Q' is given by the sum 



from the parentheses in (2.7). We call the resulting reaction system (irw(Af), Kw) 



a projection of Af onto W. Note that X is not unique. Which variables Cj, i G CVF 



are written in terms of t and which are written in terms of Cj, i G W is a matter of 
context. For instance, Example [H] below describes two different functions K associated 
with system (2.5| projected onto {1,2}. 

A natural way of defining K for K-variable mass-action systems is to include 
Ci, i G CW in the rate-constant function: K(t, c) = ~K(t)-Kw(c) 7Tw ^ P \ where K(t) = 
n(t)Tr£ W (c(t)) 7TCw ( p \ The differential equations (2.7) in this case are 



~*w(c) = Yl ( E ^P'(t)n Cw (cr^ p Aw w (cf(Q'-Q). 

Q^Q'eiZw \ {p^p'eiz-. ) 

irw{P)=Q,ir w (P')=Q'} 

(2-8) 

Note that this projection has the form of K-variable mass-action, with rate-constant 



function for the reaction Q — > Q' G IZw given by the second sum in (2.8 1. However 



this rate-constant is not necessarily bounded in a compact interval of (0,oo). 



Example 5. Following (2.8), the projection of reaction system (2.5) onto W — 
{1,2} i.e., onto species A and B, can be written in the K-variable mass-action form 
(without necessarily being K-variable mass-action): 

with differential equations 

C A = Kl(t)c B - K 2 {t)c A - K 4 (t)c A C B - K 5 (t)c A + K 6 (t) 
C B = -Kl(t)c B + K2(t)cA + K 3 {t)c A - K 4 {t)c A C B 

where (recall equation \2.(ty ) Ki(t) = Ki(t)crj(t), K2(t) — K-2{t)cc{t),~K~z{t) — n 3 (t)cc(t), 
and Ke(t) — n§(t)cr>{tf. 

On the other hand, let T — {c(t) \ t > 0} be a trajectory of \2.5ty with initial 
condition (a, (3, 7, rf) G K>o- Then c A + Cd = ol + r), c B + cc = + 7 and therefore 



7T{i j 2}(7 1 ) is a trajectory of the following projection of (2.5): 



B^A^A + B^O 2A^A^0, 

K 2 

where, denoting iri ^2} (c) = (x,y), the rate function K (t, (x,y)) is given by K\(t, (x,y)) 
Ki(t)y(a + t]-x), K 2 (t,(x,y)) = K 2 (t)x(f3 + j-y),K 3 (t,(x,y)) = K 3 (t)x(f3+j-y), 
K 4 (t, (x,y)) = Ki{t)xy, K 5 (t, {x,y)) = n 5 (t)x 2 and K 6 (t, (x,y)) = n 6 (t)(a + r\ - x) 2 . 

Remark 2.1. Let W C {1,. . . ,n] and P,P' G C. If n w (P) ^ n w (P'), then any 
directed path in 1Z from P to P' projects onto a directed path from ttw(P) to ttw(P') 
in TZw ■ If, on the other hand, ttw(P) — ^w(P'), then a directed path from P to P' 
either projects onto a cycle in TZw or * s eliminated by the projection. 

Therefore projection preserves weak reversibility: if Af is weakly reversible, then 
so is 7r(A/"). This result appears in 122}/ and is also the object of Lemma 3.4- in 

2.5. Complex-balanced systems and deficiency of a network. Complex- 
balanced systems are defined in the context of mass-action kinetics, i.e., the rate- 
constants are fixed positive numbers. 

Definition 2.10. An equilibrium G R> of a mass-action system (TZ,S,C, k) 
is called complex-balanced equilibrium if, at c*, for any complex Pq G C, the flow into 
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Pq is equal to the flow out of Pq. More precisely, for each Pq £ C we have 

P-vPo Po^P 

A complex-balanced system is a mass-action system that admits a strictly positive 
complex-balanced equilibrium. 

Definition 2.11. Let (S,C,1Z) be a reaction network withm complexes, I linkage 
classes and whose stoichiometric subspace has dimension s. The deficiency of the 
reaction network 1Z is m — I — s. 

The deficiency of a reaction network is always non- negative |12j . It has been 
shown that weakly reversible systems whose deficiency is equal to zero are complex- 
balanced [12] . This remarkable fact reveals a large class of mass-action systems which 
are complex-balanced regardless of the choice of their rate constants. 

2.6. Persistence and the sub-tangentiality condition. 

Definition 2.12. A trajectory T(c ) = {(x±(t), . . . , x n (t)) \ t > 0} with positive 
initial condition Cq G R™ of an n-dimensional dynamical system is called persistent 
if 

liminf Xi{t) > for all i € {1, . . . , n}. 

t— >oo 



Some authors call a trajectory that satisfies the condition in Definition |2.12| 
strongly persistent [2 9) . In their work, persistence requires only that lim sup^^ Xi (t) > 
for all i € {1, . . . , n}. We say that a dynamical system (or a reaction system) is per- 
sistent if all its trajectories with strictly positive initial condition are persistent. 

Definition 2.13. Let T(c ) = {x(i) | t > 0} denote a forward trajectory of a 
dynamical system with initial condition Cq € M> - The w-nm it se t of T{cq) is 

lim^ T(co) = {I € K™ | lim n _ i . 00 c(t n ) — I for some sequence t n oo}. 

The elements o/lim u T(co) are called lu -limit points ofT(co). 

Note that a bounded trajectory of a dynamical system with positive initial con- 
dition is persistent iff it has no w-limit points on cM™ . 

In this paper we will prove persistence of trajectories T(co) for reaction systems 
with special properties. Our approach will consist of showing that a certain convex 
polyhedron included in M™ contains T{cq). To this end we will use the following 
version of a result of Nagumo [8]. Recall that for a closed, convex set K C W 1 and 
for x £ K the normal cone of K at x is defined as follows: 

N K (x) = {n G E" | n • (y - x) < for all y e K}. 

Theorem 2.1 (Nagumo, 0). Let K C W 1 be a closed, convex set. Assume 
that the system x(t) = f(t,x(t)) has unique solution for any initial value, and let 
T(co) = {x(t) | t > 0, x(0) — c } be a forward trajectory of this system with c g K. 
Lf for any t Q >0 such that x(to) € dK we have the sub-tangentiality condition 

n • /(t , x(t )) < for all n e N K (x(t )) 



then T(cq) C K. 
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Fig. 3.1. (a) Essential supports corresponding to vi and V2; (b) Example of a set of inward 
vectors; (c) Illustration for the proof of Lemma \3. 6\ 

3. Lower-endotactic networks. Let Af = (S, C, TV) be a reaction network with 
species X±,... ,X n and let S C K™ denote its stoichiometric subspace. By a useful 
abuse of notation, we view the source complexes of Af as lattice points in Z™ : 

SC(Af) = {(mi, . . . , m n ) € Z> | m\X\ + . . . + m n X n G C is a source complex}. 

In this section we revisit the notion of lower-endotactic network, first introduced 
in |10j for the case of two-species networks, and we extend it to planar reaction 
networks, defined below. We let aff (Af) denote the affinc hull of C, i.e. the minimal 
affine subspace of W 1 that contains C. 

Definition 3.1. The reaction network Af is called planar j/dim(aff(jV)) < 2. 

Let aff + (AT) = aff (Af) l~l K> . The following definition is similar to the one in [lOl 
section 4]. 

Definition 3.2. Let Af — (S,C,1Z) be a reaction network such that aff (Af) has 
dimension two and let v be a vector in S. 

(i) The v-essential subnetwork Af v — (S,C V ,1Z V ) of (S,C,1Z) is defined by the 
reactions oflZ whose reaction vectors are not orthogonal to v: 

K v = {P -> P' e K | (P' - P) ■ v ^ 0}; 

C v is defined as the set of complexes appearing in reactions oflZ v . 

(ii) The v-essential support of Af is the supporting line L of comr(SC(Af v )) that 
is orthogonal to v and such that the positive direction of v lies on the same side of 
L as SC(Af v ); (in other words, for any P 6 aff (Af), the intersection of the half-line 
{P+tv\, t > 0} with the half-plane bounded by L that contains SC(Af v ) is unbounded.) 
The line L is denoted by esupp v (Af). 

Figure |3.1[ a) illustrates the notion of v-essential support for a planar reaction 
network with six complexes and four reactions. This reaction network has two source 
complexes and note that Af Vl is equal to Af, whereas Af V2 is strictly smaller than Af 
and contains only one source complex. 

We denote by esupp v (Af)<o the intersection of a.ff+(Af) with the open half-plane 
in aff (Af) bounded by esupp v (Af) that does not contain the positive direction of v : 

esupp v (Af) <0 = {P e aff + (AT) | (P - Q) ■ v < for all Q € esupp v (Af)} 

and we define esupp v (Af)>o similarly. 

Definition 3.3. Let C C K.™ be a closed and convex set, and let S be the linear 
subspace ofM. n such that the affine hull of C is a translation of S. Then 

iv(C) = - (J (JV c (x) n S) 
ii 



is called the set of inward vectors of C. Here dC denotes the relative boundary of C. 

An example of a set of inward vectors for a two-dimensional set C is depicted in 
Figure pntb). 

Remark 3.1. (i) iv(C) is a convex cone and if C is bounded then iv(C) = S. 

(ii) If C is a half-line then iv(C) consists of all vectors parallel with C pointing 
in the unbounded direction of C. If C is a bounded line segment, then iv(C) consists 
of all vectors parallel with C. 

(Hi) If aff(C) has dimension two, then the set of inward vectors iv(C) is two- 
dimensional and consists of the normal vectors v G S of all supporting lines L of C 
such that the positive direction of v is on the same side of L as C ( see Figure \3. l\ b ) 
for an example). 

Definition 3.4. Let A! be a planar reaction network with stoichiometric subspace 
S. Then Af is called lower-endotactic if the set 

{P^ P' | P G esupp v (Af) and P' G esupp v (Af) <0 } (3.1) 

is empty for all nonzero vectors v G iv(aff + (A/")). 

Definition |3.4[ ii) is easily explained by the "parallel sweep test" [10] . A reaction 
network Af — (S,C,1Z) is lower-endotactic if and only if it passes the following test 
for any nonzero inward vector v of &ff + (Af): sweep the plane aS(Af) with a line L 
orthogonal to v, coming from infinity and going in the direction of v, and stop when 
L encounters a source complex corresponding to a reaction which is not parallel to 
L. Now check that no reaction with source on L points towards the swept region. If 
TZ V = 0, then all reaction vectors of TZ are perpendicular to v and L never stops in 
the parallel sweep test. In this case we still say that the network has passed the test 
for v. 

A reaction network is lower-endotactic if its reactions with sources that are "clos- 
est" to the boundary of &ff + (Af) point "inside" aff + (Af). Note that this special 
property of lower-endotactic networks in the lattice space Z" is analogous with the 
behavior of persistent trajectories in the phase space K™ : once a persistent trajec- 
tory gets "close enough" to the relative boundary of its stoichiometric compatibility 
class S'(co), it is pushed back "inside". In this sense, the requirement that a reaction 
network be lower-endotactic appears very naturally in the context of persistence of a 
corresponding reaction system. 

Remark 3.2. Following \1Q§ . a planar reaction network Af is called endotactic 
if the parallel sweep test holds for all nonzero vectors v£S, An endotactic network 
is also lower-endotactic; the two notions coincide if aff + (Af) is bounded. 

Remark 3.3. The definition of endotactic networks has been extended in JJUj for 
networks that are not necesarilly planar, using the parallel sweep test with hyperplanes 
instead of lines (\10[ Remark 4-1])- Definition 3.4 is in fact a special case of the 



following more general definition of lower-endotactic networks: 

Definition 3.5. A reaction network (not necessarily planar) with n species is 
called lower-endotactic if it passes the parallel sweep test for any inward vector of the 
non-negative orthant M> . 

Whereas the definition above is easier to state, the more technical Definition \3.4\ 
is better suited for planar networks in the context of this paper. 

Remark 3.4. A weakly reversible reaction network Af is always endotactic, and 
in particular, lower-endotactic. Indeed, if P G esupp v (Af) for some vector v G S and 
P — > P' is a reaction of Af then P' G esupp v (Af)>o, for otherwise the fact that P' is 
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also a source complex would contradict P £ esupp v (Af) . 

Remark 3.5. If &W(Af) is one- dimensional we letV be a two-dimensional affine 
subspace ofMJ 1 such that aff(A/") C V . The parallel sweep test for Af with vectors of 
iv(V nR™ ) provides the same result as the "true" parallel sweep test with vectors of 
iv(aff+(yV*)). We may pretend that && + (Af) coincides with the two-dimensional set V 
and therefore, from this point of view, lower-endotactic planar reaction networks with 
one- dimensional stoichiometric subspace do not need a special discussion. In what 
follows, unless stated otherwise, we will assume that dim(aff (Af)) = 2. 

Note that, if aff(7V") has dimension one, then aff+(A/") contains vectors with at 
most two possible positive directions (see Remark 3.1). The following lemma shows 
that, even if aff(A") has dimension two, the parallel sweep test only needs to be 
performed for a finite set of directions v (see also [THl Proposition 4.1]): 

Lemma 3.6. Let Af be a planar reaction network with dim(aff (TV)) = 2. 

(i) If aff+(7V) is bounded, then Af is lower-endotactic if and only if it passes 
the parallel sweep test for vectors v that are orthogonal to a side of the polygon 
conv(SC(A0)- 

(ii) If &ff+(Af) is unbounded, then Af is lower-endotactic if and only if it passes 
the parallel sweep test for vectors v that are either orthogonal to a side of the polygon 
conv(SC(Af)), or are generators of the cone iv(aff + (7V)). 

Proof. The sides of conv(SC(Af)) whose inward normal vectors are in iv(aff+(7V)) 
form a polygonal line C. As in Figure [3~T| c), if aff+(A/") is unbounded, we augment 
£ with half-lines of directions given by the generators of iv(&ff + (Af)). If a vector 
v G iv(aff+(A/")) does not correspond to (i) or (ii) in the statement of the lemma, 
then esupp v (Af) contains exactly one vertex P of C. Let P P' G TZ. Since the 
inward normal vectors of the two sides of C adjacent to P belong to cases (i) or (ii) 
from the statement of the lemma, it follows that P' lies in the interior or on the sides 
of the angle ZP of £, and therefore in P E esupp v (Af)yo- In conclusion, the parallel 
sweep test holds for all v £ iv(afF+(7V)) and Af is lower-endotactic. □ 

Example 6. A few examples are illustrated in Figure 3.2 The source complexes 
are depicted using solid dots and the various lines represent the final positions of the 
sweeping lines from Lemma 3.6 Note that the reaction networks in (a) and (b) look 
the same, but, since aff+(7V) is unbounded in (a) and bounded in (b), the reaction 
network in (a) is lower-endotactic, whereas the reaction network in (b) is not. The 
same thing happens for (e) and (d) . 

Affine transformations of reaction networks. An important observation that is 
used often throughout this paper is that projections of lower-endotatic networks are 
also lower-endotactic. We prove this fact in the larger context of affine transformations 
of reaction networks. Let Af = (S,C,TZ) be a planar reaction network and consider 
an affine transformation U : aS(Af) -> M d such that U(aS + (Af)) C R> . Similarly to 
the definition of a projected network, we consider the "generalized" reaction network 
[/(AO, with reactions U(TZ) = {U(P) -> U(P') | P -> P' e 11} and complexes 
in the set U(C) which are allowed to have nonnegative real coordinates. We have 
&E + (U(Af)) = a,S(U(Af)) n M> and we may ask whether U(Af) is lower-endotactic. 

Proposition 3.1. Let Af be a planar reaction network and let U : aS(Af) — > 
WL d be an affine transformation such that U(aff+(Af)) C K> . Then, if Af is lower- 
endotactic, the planar reaction network U(Af) is also lower-endotactic. Moreover, if 
Af is endotactic, then U(Af) is also endotactic. 

Proof. We show the lower-endotactic case; the proof for the endotactic case is 
similar. Also, we assume that U has rank two; a simpler version of the argument below 
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Fig. 3.2. Examples of lower-endotactic networks - (a), (c), (e) and non lower-endotactic 
networks - (b), (d), (J). 

works if U has rank one. Let w G iv(aff + (Z7(A/"))), and let L = esupp™ (U (AT)) . Since 
U takes parallel lines to parallel lines, there exists a vector v such that U~ l (L~) = 
esupp v (Af) and 

esupp"(U{N)) <0 = U(esupp v (Af)) <0 ). 

Because U(&tt + (JV)) C aff + ([/(A0), we have iv(aff + (C7(^V))) C iv([/(aff+(A0)) and 
therefore w G iv(C7(aff+(A0). It follows that v G iv(aff+(7V)). Then, if U{P) -> 
£7(P') G £7(ft) such that t/(P) G esupp w (U{Af)) and {J(P') G esupp w (t/(7V)) <0 then 
P ->• P' G 7?., P G esupp v (AT) and P' G esupp v (H)<o, contradicting the fact that TV 
is lower-endotactic. □ 

4. 2D-reduced mass-action systems. A key ingredient in the proof of our 
main persistence result consists of studying projections of trajectories of K-variable 
mass-action systems onto well-chosen two-dimensional subspaces of W 1 . These spe- 
cial projected trajectories obey a specific type of dynamics which we call 2D-reduced 
mass-action. In this section we show that bounded forward trajectories of such dy- 
namical systems are persistent. To this end we will extend significantly the ideas from 
[TO] , where they were introduced in the context of two-species k- variable mass-action 
systems. 

4.1. Definition and comparison of reaction rates. Fix an integer n > 2 and 
let I, k be two fixed elements of {1, . . . , n} such that I < k. Let pi, qi be nonnegative 
rational numbers for i G {1, . . . , n} such that, for any i, not both pi and qi are zero, 
and such that pi = q^ = 1 and p^. = qi = 0. Denote 



Pi 
?i 



VI- 
91- 



Pi+i 

qi+i 



Pk- 



Pk+l 
Qk+i 



Pn 
In 



(4.1) 



Definition 4.1. Let * be a matrix of the form (4-D- 
(i) Let Af = (S,C,7V) be a reaction network with two species, let k 
be a piecewise differ entiable function and let a G ~ 
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>0 —r J& >0 

For all reactions P — > P' G 1Z, 



we define Kp^pi : R> X R> -> R> , 



*fp_>P/(*,x) = K p^ P ,(t)(^f p+a . (4.2) 

The reaction system (Af, K) is called a 2D-reduced planar mass-action system and is 
denoted by (Af, \E I , k, a). 

(ii) For each s G {1, . . . ,p}, let Af s — (S,C s ,7t s ) be a two-species reaction net- 
work and let (Afs,^, K s ,a s ) be a 2D-reduced mass-action system. The sum (Af,K) 



Us =1 (A/" s , ^, K a , a s ) (recall Definition 2.1) is called a 2D-reduced mass-action system. 

Therefore the concentration vector c(t) — (x(t),y(t)) of a 2D-reduced mass- action 
system [j p s=1 (Af s ,^, 

n s ,a s ) satisfies the following differential equation: 

^) = E E K S , P ^p'(t)(Mt)f p+a '(p'-py (4.3) 

s=i p-»P'eR,, 

Note that, by definition, k needs not be bounded away from zero and infinity, as is 
the case for At-variable mass-action systems. However, we will require this condition 
to prove persistence of 2D-reduced mass-action systems in Corollary |4.1| 

The goal of this section is to study the persistence of 2D-reduced mass-action 
systems. One important component of our analysis is highlighting the reaction whose 



rate at time t > "dominates" the other reaction rates. In view of (4.2) we then 
consider, for A > and for any s € {1, . . . ,p], the sign of the difference 

(*x)* p+a = - A(*x)* p ' +a % (4.4) 

for all pairs of distinct source complexes P, P' of Af s . For simplicity, and without loss 



of generality, we assume that k = 2 and 1 = 1, Then (4.4) has the same sign as the 
following expression, which we denote by p(x,y) : 

A^p(x, y) = x a y-?( P3 x + q 3 y)^ a -^ . . . (p n x + q n y) p " a - q " - A D . (4.5) 
Here (a, — j3) — D(P — P') and D denotes the least common denominator of all 



nonzero pi and qi, i € {3, . . . , n}. Note that all the exponents in (4.5 ) are integers. 

The geometry of the curves a(x,y) = within IR^q is very relevant to our 
discussion. An immediate goal, which we pursue next, is to find simple approximations 
for these curves. We will see that within appropriate subsets of M> i b( x 'V) = 
may be approximated by power curves y = C'x T that are ordered in a useful way, as 
we will explain later in the paper. Let 

A= U + E^J a ' ^= (^ 1 + E*^ and A = A -5. (4.6) 
Lemma 4.2. Suppose A^O. 

(i) If ' aj3 < then for all x > there exists a unique y > such that A^ p(x, y) = 

0. 

(ii) If a/3 > then for all small enough x > there exists a unique y > such 
that A£ i/3 (x,y) =0. 

Proof, (i) If af3 < 0, without loss of generality we may take a > and /3 < 0. 
For any fixed x > 0, A^ ^(a;, y) is a polynomial in y whose coefficients are all positive, 
except for its free term —A D . The Descartes rule of signs implies that this polynomial 
has a unique positive root. 
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(ii) If a/3 > 0, we may assume that a > and /3 > 0. Jx, y) — implies 
A D yf } (p 3 x + q 3 y)^ . . . ( Pn x + q n y)^ = x a (p 3 x + q 3 yf* a . . . (p n x + q n yf^ . (4.7) 

We rewrite this equality by excluding the factors of zero power and merging the powers 
of x in the left hand side and the powers of y in right hand side. We denote 

A K = 1 + Pi and x y = 1 + ( 4 - 8 ) 

i£{3,...,n\ ££{3,...,rt| 
qi=0} Pi =0} 

and we let £/, I G {1,...,/} be the indices for which both p^ and q^ are strictly 
positive. Then we have 

A'y x y^p n x+q n y)^ . . . ( PlI x+q tI y)^0- x ^( PtiX+qii y)P^ . . . (p^x+q^y)^ = 

(4-9) 

where A' — A D Oie{i n}\{ii ii} iT^Pi We denote the polynomial in (4.9) by 
F(x,y). _ ' 

If A y /3 > A — X x a then equation (4.9) yields 

F(x,y) = d Q y K +(d 1 x)y S ~ 1 + ... + @ T _ Xy0 x^-^)y^ 

where the coefficients Cfe and Ck are positive and are obtained from expanding the 
first, respectively second term of the difference (4.9), and if \ y f3 < A — X x a we have 

F(x,y) = doy^+id, x)y^ + ... + (^_ (S ^ Q+1) ^-(A-A, Q+ i) )y S-A, Q+1 + 
+ (d=_ Xy0 x^-C^ x *-^)y^ 

- (Ca-^p-i) x^-^-^y^- 1 ... (C^ x^)y -C^x K = 0. 

In the first case, for any fixed x > 0, the coefficients of F x (y) — F(x, y) viewed as 
a polynomial in y change sign exactly once. It follows from the Descartes rule of signs 
that F x (y) = has a unique positive solution. In the second case, the coefficients 

of F x (y) that are binomials in x are of the form Cx A ~ k — Cx k and for small x 
are all either positive if A < A or negative if A > A. Therefore for small enough x 
the polynomial F x (y) changes the sign of its coefficients only once either at y A ~ x * a if 
A < or at y A «' 3 ~ 1 if A > 0. It follows that the equation F y {x) = a unique positive 
solution. □ 

implies that for a/3 < 0, the curve {(x,y) G K> | 



Remark 4.1. Lemma 



4.2 



A^p(x,y) = 0} is the graph of a function y&p ■ K>o ~ ^ IR>o- It is easy to see 
that lim^oVa p( x ) = 00 an< ^ um x->-oo Va p( x ) = 0- O n the other hand, if a/3 > 0, 
the function y^p is defined only for small x > 0: there exists ^ > such that 
{(x,y) G (Q,M£ p ) x R >0 | A* p {x,y) = 0} is the graph of y£ p : (o,M* p ) -> K >0 . 
We claim that in this case we have \im x ^,o y^s — 0. Indeed, suppose A < 0. If for 
some < I < oo , (0,1) is a limit point of the curve {(x,y) G R>o I b{ x iV) = 0}' 
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then plugging (0, 1) into yields I = 0. It remains to check that (0, oo) is not a 

limit point of the curve above. If {(x n , y n )}n>o is a sequence of points with pos itive 
coordinates such that ]mx n -^ 00 (x n ,y v ) = (0, oo) and g(x n ,y n ) = 0, from (4-9) we 
get 0{x^ a ) = 0(y~ A+x * a ) as n -> oo, which contradicts A < 0. The case A > 
follows from symmetry. 

Lemma 4.3. Suppose a/3 > 0. 

(i) If A = £/iere exist positive constants 7„ 1: ■ ■ ■ , 7^ at /or some integer N > 1 
smc/i £/iat 



a,/3 



(x, y) = /or some (x, y) e K.> i/ and only if y = 7„ ^x /or some i € {1, . . . , N}. 



(ii) If A ^ i/iere exists a strictly incr easing function r 
for function y& g introduced in Remark 4-1 



the limit 



lim 



exists, is positive and finite. We denote this limit by C^ s . 



Proof, (i). Dividing (4.9) by x A = x A and letting 7 = y/x yields 



A'^'Pfa + qill )^ . . . [p h + q tll Y^ - ( Ptl + g il7 r ia • ■ • (P« + <Hrl) Pi ' a = 0- 



We denote by F(-j) the polynomial above. The positive term in the expression of F 
has degree A, and the negative term has degree A — X x a, therefore lim 7 _ i . 00 F(j) = 00. 
Since F(Q) < 0, there exists at least one positive root of F. Denoting the positive roots 
of F by "fall ■■■ )7a iv completes the proof. 
(ii) We know from see Remark 



4.1 



that (0,0) is a limit point of {(x, y) g 



I 2 



A^ g (x,y) — 0}. Lemma 4.2 implies that this curve has a unique Puiseux expansion 



in a neighborhood of (0, 0) (see [2"5] for a discussion of Puiseux expansions). By 
making M£ 8 from Remark 4.1 as small as necessary for the Puiseux expansion to 



hold in x G (0, M* p ), we have, for all x € (0, M^ g ): 



\X 



higher order terms in x 



(4.10) 



where C£ p € K and T > is a rational number. The exponent T in rt4.10| is equal 



to the negative of one of the slopes in the lower boundary of the Newton polygon of 
the polynomial F(x;y) defined in (4.9); (see [35] for more details). Since F is the 
difference of two homogeneous polynomials of degrees A and A, its Newton polygon 



can be easily illustrated (see Figure 4.1) and the slopes of its lower boundary are 



-1,- 



A - (A - \ y l3) 

Xy/3 



-L-E 



A — (A — X x a) 



(4.11) 



if A > o r A < 0, respectively. Note that T cannot be equal to one, for otherwise 
from (4.9) it would follow that 0(x A ) = 0(x A ), contradicting A ^ A; therefore -T 
is given by the fractions in (4.11 ). Then, if we define the function r : M>o — > K>o as 
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(0,A) 




(Xy/3,A- Xy/3) 



(0,A) 




(A,0) 



(A,0) 



Fig. 4.1. Monomials of F(x; y) (see \4..{fy . The lattice point (a, b) represents the monomial x b y a 
and the lower boundary of negative slope of the associated Newton polygon is the thick polygonal 
line. The picture on the left corresponds to A > and the one on the right corresponds to A < 0. 



follows: 



t(ct) = I 1 



A I, 



^ n n 



i=l i=l 



if cr < ^qi/^Pi, 

i=l i=l 
n n 

if cr = ^Qi/^Pi, 

i=l i=l 
n n 

if cr > ^Qi/^Pi, 



i=l 



i=l 



we can easily see that T = r(a//3). Note that t is continuous and strictly increasing. 
Moreover, the value of a for which r(cr) = 1 corresponds to A = 0, and thus the 
statement in part (z) of the lemma is incorporated in part (ii). □ 

4.2. The domination lemma. The following key result reinforces our motiva- 



tion for considering differences (4.4) of reaction rates. Roughly speaking, it shows 



that if, at time t, the rate of a reaction Pq — > Pa dominates all the other reaction 
rates, then the reaction vector Pq — P "dictates" the direction of the flow c(t). 

Lemma 4.4. Let {Af,K) = (S,C,1Z,K) be a reaction system, let P — > Pq £ 
TZ and let v be a vector such that (P^ - P ) ■ v > 0. Also let U C SC(Af)\{P }. 
There exists a positive constant /i such that if for some (t, x) £ K>o x R> we have 
Kp ^ P i (i, x) > and 



K Pa ^pi t (t, x) > ptKp^p, (t, x) for allP^ P' £K with P e U 



then 



(Kp ( ^p i (t,^(P[ ) -Po)+ ]T K P ^ P ,(t,x)(P'-P)) v 



> 0. 
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Proof. We take 



|vii w p - p '\ 



A* = 



and we have 

(A> fW p Q (i,x)(^-P )+ £ W(t,x)(P'-P) -v 

>K Po ^ Pi (t,x)(P' -P )-v- ((V/A^-^fr*) l|v||||P-P'|| ) 

p^P'.Peu \ ) 

> ((P^-Po)-v-(l//i)||v|| £ ||P'-P||W^(f,x) = 0, 
V p^P'en ) 

where the first inequality was obtained using the Cauchy-Schwarz inequality. □ 

Remark 4.2. If (S,C,lZ,ty,n,a) is a 2D-reduced planar mass-action system 
and if n(t) € (r/,l/r/) K for some t > then K P ^ P ,(t,x.) = K P ^ P ,(i)(1'x)* p+Q , 
and therefore (*x)* p °+ a > (^/?/ 2 )(*x)* p+a implies K Po ^ P >(t,x.) > ^iK P ^ P ,(t,x.) 



which is exactly the condition needed in Lemma 4-4 Or, using the notation in (4-5) 
and letting (a, —j3) = D(Pq — P) : 

A a$ 2 ( x ) > lm P Ues Kpo^M^) > ftKp^ P >(it,x). (4.12) 



4.3. Geometric constructions in the phase plane. Our strategy for proving 
that a trajectory T(cq) of a certain reaction system is persistent relies on building a 
convex set C + C R> that contains T(co) and stays away from cM> . As in [TP] , 
we partition the phase plane into subsets where one reaction rate dominates all the 
others, and therefore, by Lemma |4.4| its corresponding reaction vector dictates the 
direction of the vector field. The set C + is constructed such that, on each subset 



of the partition, the dominating reaction vector (and therefore, by Lemma 4.4 the 
vector field) , points towards the interior of £ + . This is the rather simple idea behind 
the proof of Theorem |4.1[ but the technical details involved are quite delicate. We 
start with the construction of the set L + , which is discussed next. 

For any s € {l,...,p}, let N s = (S,C S , TZ S ) be a lower-endotactic two-species 
reaction network and let (N : K) — {J^ =1 (Af s , K s ,^,a s ) be a 2D-reduced mass-action 
system. We also let r\ < 1 be a positive constant. Let D denote the least common 
denominator of all nonzero elements of vp. For each s € {1, . . . ,f>} let 

{rf, . . . , r e s (s) } = {a//3 \ (a, -0) G {D(P - P') \P,P'e SC(M s ),af3 > 0}}; 

we assume that rf < . . . < r*^ a j and define the set 

V = {(l,f 1 ),...,(l,f e )}=U{(l,^)|ie{l ! ... ! e( S )}} (4.13) 

where we take fi < . . . < f e . Let {i, j} be the standard basis of the cartesian plane. 
Since Af s is endotactic, for each s € {1, . . . ,p} and for all vectors n £ VU {i, j}, there 
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exists a reaction 

P s ,n -> P' s „ € ^ a such that P Si „ g esupp n (Af s ) and P s ! n g esupp n (Af s ) >0 . (4.14) 



Note that there might exist multiple reactions as in (4.14), out of which P s 
P' s n is chosen and fixed for the remaining of this paper. 

If /i^,, denotes the constant from Lemma 4.4 that corresponds to reaction P s 



P s ' n v = n and U = SC(Af s ) D esupp n (J\f s ) >0 we define 
fi = max {/i„ jS | n g V U {i, j}, s g {1 



(4.15) 



We also let 



V = Q {D(P - P') | P, P' e SC(Af s )}, 



s = l 



and finally, inspired by Remark 



4.2 



we denote A = n/rj 2 . 
Let M > 1 be a fixed number. We choose < 8 < 1 and < £ < 1 to satisfy the 
following properties: 

1 



(PI) 



{C£/>,7£< I (a,-/3)g2W>0}, 
> | max{C^, 7 ^ | (a, -/?) g 23, a/3 > 0} 



5 < — mm 
1 
8 



4.3 



where C A a and "f A 4 are defined in Lemma 
(P2) all pairwise intersections from the strictly positive quadrant of the 2e curves 
y = Sx T ^\ y = (l/8)x T ^\ ie{l,...,e}, lie in (£,oo) 2 ; 

(P3) the square [0,£] 2 lies below the curves A A g(x,y) = for all (a,—/?) g 
V with a/3 < 0; 

(P4) for all (a, -0) g V such that a,3 > we have £ < (recall that M A fi 

is such that : (0,M£g) admits a Puiseux series representation), and for all 
2; g (0, £] we have 



Clearly, £ can be chosen small enough so that (P2) and (P3) are satisfied (see 



Figure 4.2 ). The existence of £ that also satisfies (P4) is a consequence of Lemma 
Since r is a strictly increasing function, condition (P2) implies that for x £ (0, £ 
have 

(l/8)x T ^+^ <Sx T ^ for all i g {1, . . . , e - 1}. 
Moreover, for x < £, (PI) and (P4) imply that for a/3 > the curve {{x,y) g 



4.3 



We 



4.2|. 



A A g(x,y) — 0} lies between the curves 8x T ^ a ^^ and (l/<5)x r ( a /^ (see Figure 

Now, for the actual construction of C + , first assume that V ^ and start with 
a point A g {0} x (0, £] on the j/ axis. We choose a point A 2 between the curves 
y = [\/8)x tIj ' 2 ^ and y = 8x T ^ ri ^ such that the slope of the line AqA 2 is — 1/fj.. 
Inductively, choose points Aj+i, i G {2, . . . , e — 1} such that 

{%) the point Ai + \ lies between the curves y = [1/ 8)x Tt ^ ri+1 ^ and y = 8x Tlj ' i ^ 
(it) the line AiAi + i has slope — 1/fj. 
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Fig. 4.2. Construction of C + . Each dashed line represents a curve g = in the positive 

quadrant and is labeled with a relation that the corresponding a and /3 satisfy. The curves y = 5x T ( r ' ) 
and y = (l/S)x T ^ ri ^ are depicted using solid lines. The highlighted area represents the set £ + . 

Finally, A e+ 2 is defined on the x axis such that A e A e+ 2 has slope —l/f e . 

The polygonal line [^4o^2^3 ■ ■ ■ A e ^iA e A e+ 2\ is convex. We move Aq closer to the 
origin if necessary, such that all the points Ai defined above lie in the square (0,£] 2 . 

If V = we let Aq and A3 be two points on the y and x axes, respectively, such 
that Ao, A^ € [0,£] 2 and the slope of the line A0A3 is -1. 

The last step of the construction consists of defining d small enough such that 
(P5) and (P6) below are satisfied: 
(P5) If V ^ then the vertical half-line {d} x [0, 00) intersects the segment Aq A2 at 
a point A\ above the curve y = (1 / '8)x T ^ r ^ '; also, the horizontal half-line [0, 00) x {d} 
intersects the segment A e A e+2 at a point A e+ i below the curve y — 8x T ( r "K If 
V = then d is chosen such that d < £/2, and the intersection of the vertical half- 
line {d} x [0,oo) with the segment j4o^3 is denoted Ai, whereas the intersection of 
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the horizontal half-line [0, oo) x {d} with the segment ^0^3 is denoted 
(P6) A£ iP (d,M) > for (a,-f3) € V,a < 0,/3 > 0, 

A a!/3( M - d ) > for («. -/?) e a > 0, /3 < 0; 
Recall that M > 1 was fixed at the beginning of the construction. It is easy to 



see that (4.5) implies that (P6) holds for d small enough. 
Let 

C = {A 1 + (Q,t) |t>0}uL4!...yl e+1 ]u{A e+1 + (t,0) | t > 0} 

be made out of the polygonal line \A\ . . . A e ], (called, for future reference, the finite 
part of C) completed with a vertical and a horizontal half-lines (the union of which 
we call the infinite part of £). Finally, we define C + = conv(£). 

To indicate the quantities that C + depends on, we write C + = C + ({Af s }i< s < P , Tj, M). 
The polygonal line L4 , . . . , A e+2 ] will also be useful later in this paper. We denote it 
by C — £({Af s }i< s <p, rj) (note that M is not required for C) and we let £ + be the 
unbounded part of R> that is delimited by C. 

4.4. Persistence of 2D-reduced mass-action systems. The following theo- 
rem is the main result of this section. 

Theorem 4.1. Let (N,K) — \J P S=1 (N S , k„, a s ) be a 2D-reduced mass-action 
system where Af s = (S,C S ,TZ 8 ) is lower- endotactic for all s £ {1, . . . ,p} and denote 

/(t,x) = J2 E k >w (t)W p+as (p' p). 

s=i p->P'eK s 

Let ij £ (0,1) and M > 1 be real numbers. Then for any t > such that K s (t) £ 
(t], 1/?7) Ks for all s G {1, . . . ,p}, and for any x G C({Af s }i< s < p , 17, M) n [0,M] 2 we 
have 

n • f(t, x) > for all n G -N c + (x). (4.16) 
Proof. The cone —N c + (x) is degenerate unless x is a vertex of C and its generators 



belong to V U {i,j}, where V is the set of vectors defined in (4.13). It then suffices 



to show that for any s G {1, . . . ,p} and for any n G —N c + (x) n (V U {i, j}) we have 
( Ks,P^P'(t,x)(P'-P))-n>0 (4.17) 



for all s G {1, . . . ,p}, where K S! p-ypi (t, x) = K s .p^p'(t)(5 , x) ,I ' p+as is the rate of the 
reaction P — > P' G 1Z S at time t. We fix s G {1, . . . ,p} and recall that lZ s , n = {P — > 
P' G 1Z S j (P' — P) ■ n ^ 0} denotes the n-essential subnetwork of 1Z S . The inequality 



(4.17) is trivially true if 7\L s , n = 0. Otherwise, recall the reaction P s n — ¥ P' s n G 1Z 



from (4.14). We rewrite the left hand side of (4.17) by separating the reactions with 
source on esupp a (J\f s ) and emphasizing the reaction P Si „ — > P' s n : 
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( K s>P ^ P ,(t,x)(P'-P)Yn=l K a> p^ P ,(t,x)(P' -P) + 



P£esupp™ (J\f s 
P-j.p'^P s „^P^ n } 



-K atP , tB ^p l Jt,x){P'. in -P., n )+ K s ^ P ,(t,x)(P' -P) n. 



{P^P'GAf s ,„| 

P^esupp n {J^ B )} 



Since all source complexes of 1Z SiI1 lie in esupp n (M s )>o, the reaction vector P' — P 
with source P e esupp n {N s ) satisfies (P' — P) • n > 0. It is therefore enough to show 
that 



K s,P„^PL a (*,x)(i* -Pn)- 



E 

{p^p'eK s ,„| 

P^esupp n (AA s )} 



X s , P ^(i,x)(P'-P) n > (4.18) 



in order to verify ( |4.17 1. In turn, (4.18) will follow from Lemma 4.4 with U 
SC{N ' ) n esupp n (lZ s ) > o and the fact that 



K P s n ^P' n (^ x ) > ^Kp^ P ,(t,x.) 



(4.19) 



for all P — >• P' e 7?. s with PeW (recall /z from ( 4.15[ )). Therefore showing inequality 
(4.19) will complete the proof of the theorem. 



As noted in Remark 4.2 (4.19) is implied by 



A£„(x) > 0, 



(4.20) 



where (a, — 0) = D(P S I1 — P) and A = fi/i] 2 . To verify (4.20), we consider different 



cases, according on the location of x = (x,y) on C. First suppose that x lies on the 
line segment [AjAj+i] for some i € {1, . . . , e}. Then (see Figure 4.2), if fy < fj < ry/, 
or equivalently, if z' < i < i" , we have 



(l/<5)x T ('>> < y < <S:r r ^. 



(4.21) 



Depending on the sign combination of a and f3, the source complex P may belong to 
one of the three shaded regions in Figu re |4~3{ a) . 



Region I. Here a < and (3 < 0; 14.201) is equivalent to y > yt n( x )- The set 



{r? | j G {1, . . . , e(s)} with t) > f<} 

is nonempty, since a//3 is one of its elements. If r| o denotes the minimum of this set, 
then r| o < a/ (3 and it follows that 

y > (l/S)x< r V > (3/2)C^x^ > y^(x). 
The first inequality is implied by (4.21), since > fy The second inequality holds 



because of (PI), because x < £ < 1 and because r is increasing. The last inequality 
is a consequence of (P4). Note that the calculation above corresponds to A ^ 0, but 
the same argument works if A = with g replaced by 7^ = max 7^. 
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Region II. We have a < and /3 > 0. From (P3) we know that x £ [0, £] 2 is below the 
curve Af j(3 = 0, and so A^ /3 (x) > 0. 
Region III. This case is similar to Region /. 

Finally, suppose that x lies on one of the two unbounded sides of C, for instance, 
on the vertical side. Then a < and there are two regions for P, according to the 
sign of (3 (see Figure |4~3{ [b)). For region / (/3 < 0), the proof is the same as in the 
case of region / from Figure |4~3^ a). For region II we have j3 > and a < and, since 
y < M, we have 

A^(x)=A^(d )y )>A^(d,M)>0 

from (P6). □ 

The following corollary follows from Theorem |4.1| using Nagumo's Theorem |2.1| 
and implies that bounded trajectories of 2D-reduced mass-action systems are persis- 
tent. 

Corollary 4.1. Let Ug =1 (A' r s ,K s , v I',a s ) be a 2D-reduced mass-action system 
where Af s — (S,C S ,TZ S ) are lower- endotactic networks for all s € {1, ...,p}, and 
suppose that K s (t) £ (r],l/ri)' R ' s for all s £ {1, . ..,£>} and all t > 0. Let T(cq) 
be a trajectory of (J\f,n) and let M > 1 be such that T(cq) C [0, M] 2 . If C + = 
C + ({Af s }i< s <p, ^P, rj, M) is constructed such that cq £ C + , then T(cq) C C + . 



Remark 4.3. Corollary J^.l remains valid if instead of cq £ £ we have cq £ 
C = C {{Af s }i<s<p, ^, ?])■ In that case we conclude that T(co) C C . 

We conclude this section with the following result, which will be useful in section 

El 

Corollary 4.2. Let [J p s=1 (J\f s , K Sl ^/,a s ) be a 2D-reduced mass-action system 
where J\f s — (S ,C Sl TZ s ) are lower-endotactic networks for all s £ {1, . . . ,p} and sup- 
pose that K s (t) £ {r],l/rf) R ' a for all s £ {l,...,p} and all t > 0. Then there exist 
e > and r > such that if cq — (xq, yo) £ [0, e] 2 then x + y > t(xq + yo) for all 
(x,y) £ T(co). _ 

Proof. Let L = L({N s }i< s < P , *, rj) and let e > be such that [0, e] 2 n£({7V s }i< s < p , 
0. Once C is constructed, we can shift it as close to the origin as desired. In partic- 
ular, if cq £ [ 0, e] 2 , we ma y con struct Cq — £o({A/ r s }i< s < p , 4 r , rj) such that cq £ Co 



!-' — 



Corollary 4.1 and Remark 4.3 imply that T(cq>) lies in the unbounded part of ii^ >0 
delimited by Cq. We draw lines through c that are parallel to the extreme line seg- 
ments of Cq and denote their intersection with the coordinate axes by (0, jj) and (x, 0) 



(see Figure 4.4 1. Then T(cq) lies above the line x + y = min{i, jj}. Using the notation 
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from (4.13 1 we have x = xq + Tiyo and y = Xo/r e + yo and the conclusion follows by 



choosing r = min{l, n, l/r e }- □ 

5. Persistence of K-variable mass-action systems with two-dimensional 
stoichiometric subspace. The main persistence result of this paper is the following. 

Theorem 5.1. In any n-variable mass-action system with two-dimensional sto- 
ichiometric subspace and lower- endotactic stoichiometric subnetworks, all bounded 
trajectories are persistent. 

A little additional terminology and a couple of lemmas are needed to arrive at 
a proof of Theorem |5.1| For the remainder of this section, we fix a k- variable mass- 
action system (J\f,n) — (S,C,1Z, k) having n species, stoichiometric subspace S of 
dimension two, and lower-endotactic stoichiometric subnetworks Af s = (S,C S ,1Z S ), s £ 
{1, ...,p}. We assume that «(t) £ (ry, l/r/) n . We also let T(c ) = {c(f) = (xi(t), . . . , x n (t)) 
t > 0, c(0) = Co} be a bounded trajectory of (A/ - , k) such that Co G M™ and 
T(c ) C [0,A/]™ for some M > 1. 

5.1. Preliminary setup. For W C {1, . . . , n} we let 

Z w = {(xi, . . . , x n ) e R n \xi = for all i e W}. 

The relative boundary of the polyhedron S(cq) is included in 9M™ and we may 
identify a face of S(cq) by the minimal face of M™ that contains it. More precisely, if 
a face of S(cq) is included in Z\y and W C {1, . . . , n} is maximal with this property, 
then we denote that face by Fw- Note that if W 1 C W, then Fw C i^w 1 ; an d 
F w = Fvvi if and only if W = W 1 . 

Decreasing n if necessary, we can assume that 5* intersects the open positive or- 
thant, S n M™ ^ 0. Indeed, if this is not true, then some coordinates Xi, i 6 ZV C 
{1, . . . , n} of c(t) are constant and may be disregarded, by replacing (Af, n) with its 
projection (A/y, Ky) onto V. Note that the properties of N are inherited by AV : the 
stoichiometric subspace of Ay has the same dimension as S; the projected kinet- 



ics is k- variable mass-action (see (2.8)) and Ay has lower-endotactic stoichiometric 
subnetworks from Proposition 3.1 finally, 7ry(T(co)) € [0,M] y . 

Each vertex Fw of S(cq) has two adjacent edges, which we will henceforth denote 
F w i and F w %. We have W 1 U W 2 C W and also W 1 n PF 2 = 0, for otherwise 5 
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would not intersect the interior of the positive orthant. Let v 1 (W) = (pi, . . . ,p n ) and 
v 2 (W) = (qi,.-.,q n ) be vectors (unique up to positive scalar multiplication) along 
F w i and F W 2 respectively, such that the cone with vertex at Fyy generated by F w i 
and F\^2 contains S(cq). Then, for any i £ {1, . . . , n}, not both pi and qi can be zero, 
for otherwise, again, S C 5R> . 
Remark 5.1. We have 

Pl = 0fori£ W 1 ,p i ^Qfori£ ZW 1 and Pl > for i £ WXW 1 , (5.1) 
qi = 0forie W 2 , qi ^ for i £ ZW 2 and q % > for i £ W\W 2 . 

Indeed, if pi = for some i £ CW^ 1 , then F\yi — Fyyi\ju\, contradiction. On the 
other hand, for any x = (xi, . . . ,x n ) £ F\yi\Fw, we have x = Fw + a x (Pij ■ ■ • ,Pn) 
for some a x > 0, and therefore Xi — a x pi for all i £ W. Since Xi £ Fyyi it follows 
that Pi = for i £ W . Moreover, since Xi > for i £ W\W , we have Pi > for 
i £ W\W . The explanation for qi, i £ {1, . . . , n} is similar. 

Remark 5.2. If (k,l) £ W 1 x W 2 , then we may rescale ^(W) and v 2 (W) such 
that pi = (ft = 1. Moreover, by swapping v 1 (M^) with v 2 (W) if necessary, we may 
also assume that I < k. Let 

Q = ( V 1 (W) v 2 (VK)) (5.2) 

be the matrix with columns v 1 (W) and v (W). Since v 1 (W) and v 2 (W) generate S, 
we have 

x = 97Tfc i /(x) for any xeS (5-3) 
(here, j/x = (xi, . . . , x n ), then 7rfc,j(x) = [xi, xu) )■ In particular, it follows that 

x - F w = 67rfc,;(x) for all x £ S(c ). (5.4) 



Remark 5.3. Let (k,l) £ W 1 x W 2 and let N s = (S,C s ,K a ), s £ {l,...,p} 
denote the stoichiometric subnetworks of Af. Then, for each s £ {1, . . . ,p}, we may 
choose a s € E™ such thatC s C S + a s and irk,i{a s ) = (0,0) (recall thatiTk,i denotes the 
projection onto {k, 1} ). Indeed, suppose C s C S + (71, . . . , J n ) for some (71, ... , j n ) £ 
R". Then we look for a s such that 

a s -( 7l ,..., 7 „) £ S^s V sm{w 1 (W),w 2 (W)} 



and TTk,i(a s ) = (0,0). Assuming pi = = 1 as in Remark 5.2, we define 
a s = (71,.. . ,7„) - 71V 1 (W) - 7fcV 2 (W). 



Example 7. It might be helpful at this point to illustrate the notations in- 



troduced thus far in this section by revisiting the network in Example 2.1 If we 



let Co = (1,1,1,1), it is not hard to see that S(cq) is a square with vertices at 
(0,0,2,2), (0,2,0,2) (2,0,2,0) and (0,0,2,2). Let W = {3,4} and consider the 
vertex Fw — (2,2,0,0) of S[cq). For any {x\, X2, X3, 2:4) £ S(cq) we may write 



( Xl ,x 2 ,x 3 ,x 4 ) - (2, 2, 0, 0) = x 3 (0, -1, 1, 0) + x 4 (-l, 0, 0, 1) 
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and therefore V 1 (W) = (0,-1,1,0), v 2 (W) = (-1,0,0,1) and W 1 = {4}, W 2 = 
{3}. The face F w i of S(cq) is {(2,2 — x, x, 0) £ M \x £ [0,2]} and is parallel to 
v (W). As for the vectors a s from Remark 5.3 corresponding to our two stoichiometric 
subnetworks, we have a\ = (2, 2, 0, 0) and ai = (2, 0, 0, 0). 

Remark 5.4. Since v 1 (W / ) and v 2 (W) generate S, it follows that 

x = 07Tfc j ;(x) for any x £ S. (5-5) 

In particular, we have ker^; (1 S — and therefore ir^^ is injective on S. It follows 
that TTk,i is also injective on a s + S. 

If Q is a column vector in R^ k ^ ~ M 2 then <dQ £ S and we have 

7Tfc,i(a s + 9Q) = Tr k ,i(a s ) + Tr kt i(QQ) = Q. 

Since, as explained above, TTk,i is injective on a s + S, we conclude that 

^}(Q)n(a s +S) = {a s + eQ}. (5.6) 



5.2. A glimpse into the rest of section[5} Before we dive into the technical 
arguments that lead to a proof of Theorem |5.1| it is worth considering a couple of 
examples. The aim is to illustrate how the machinery of projections and 2D-reduced 
mass-action systems comes into place and to hint at the idea behind the proof of 
Theorem 15.11 

Let us first consider the following ^-variable mass-action system: 

A + B^2C A + C^A. (5.7) 

We denote x = (c^, cb,cc) the concentration vector and we write the correspond- 
ing differential equations in the form 

cfj = Kl (t) X A+B ^-lj +n 2 (t) X 2C ^ 1 j + K3 (i)x A+c ^ j + K4 (i)x A |o 

(5.8) 



Note that the stoichiometric subspace of (5.7| is S = {(21, X2, X3) £ M 3 \xi = X2} 
and it intersects the positive orthant M> - Let then cq £ Sflf^ be a positive initial 



condition for (5.8) such that T(cq) is bounded. If we let 




then for x = (ca 7 cb, cc) £ T{c$) and denoting y = (ca, cc) we have 

x = fy. (5.9) 

Moreover, note that 

a+b= 1 1 1 = *(:j =*(A), 




where, by a useful abuse of notation, the argument A of iff is viewed as the vector of 
two coordinates (1,0) in the plane AC. One can write similar equalities to obtain 



A + B = S>(A), 2C = ff(2C), A = iff(A) + (0, -1,0)*, A + C = iff(A + C) + (0,-1,0)*; 

(5.10) 

denote a\ — (0, 0, 0) and a 2 = (0, —1, 0). 

Projecting ( |5.7[ ) onto {A, C} yields the reaction network 

A^2C A + C^ A 



whose dynamics is obtained by substituting (5.9) and (5.10) into (5.8): 



y = Kl (t)(*y)*(' 4 )+^ ( f\ + K2 (t)(*y)*( 2C )+^ f ^ (5.11) 

+ K3 W(*y)* (A+c)+a2 (V) + <^(t)(*yf {A)+a2 (J) • (5.12) 

This kinetics corresponds to the 2D-reduced mass-action system (A/"i, iff, re', ai) U 
(7V 2 ,*,K",a 2 ), where 



Ni = {A ^ 2C}, JV 2 = {A + C ^ A}, 



and 



Corollary 4.1 implies that the trajectory of y is persistent; from (5.10) it follows that 



T(co) is persistent as well. 

The argument above can be written in the general case without much additional 



effort; this is done in Proposition 5.1 Note that, although it illustrates very well the 



use of projections and 2D-reduced mass-action, the example discussed above is rather 
special by insisting that S(cq) contain the origin. To see what issues might arise if this 
is not the case, let us next revisit the system (2.5), which we assume to be re- variable 



mass-action. Choose cq in the same stoichiometric compatibility class as (1, 1, 1, 1). 
Since S(cq) is bounded, so is T(cq). Note that for any x = (ca, cb, cc, cd) G T(cq) we 
have ca + c D = c B + c c = 2. 

We will now give an heuristic explanation of the fact that ca cannot become too 
small (the same reasoning may be applied to the rest of the concentrations). We 
aim, as in the previous example, to project our system onto a 2D face of IR> and 
realize the projected dynamics as 2D-reduced mass-action, in order to conclude that 
Cyi(t) stays bounded away from zero. Let us consider the projection onto {A, B}. As 
illustrated in Example [5j the projected network 



B^A^A+B^O 2A^A^0 



(5.13) 



is lower-endotactic and the projected dynamics can be written in the form 
' ' 1 «i(*)cb ( \ J +K 2 (t)c A f ) +K 3 (t)c A K 



CB 



(5.14) 



K 4 (t)c A C B 



K 5 (t)c A 



K 6 (t) 
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where 



K 1 (t)=K 1 (t){2-C A (t)), 
« 4 (t) = K 4 (t), 



K2(<) = N2(*)(2 - C fl (t)), 7Ss(t) = K 3 (t)(2 - C B (t)), 

K 5 (t) = K 5 (t), K 6 {t) = n 6 {t){2 - c A {t)f . (5.15) 



This is a 2D- reduced mass-action system (in fact, this would be k- variable mass- 
action system, pro vided we knew that Kj are bounded away from zero). Since T(cp ) is 



bounded, Theorem 



4.1 



implies that there exists a set C\ B C 



,{A,B} 



as in section 



4.3 



such that the projection of Cq onto {A, B} lies in C\ B and whenever the phase point 



(cA{to),CB(to)) of (5.14) is on the boundary of C\ B , the vector field points inside 



C\ B ; this, provided /tj(to) belongs to a certain interval away from zero and infinity. 
By inspecting the rates Ki in (5.15), we see that this condition is equivalent to saying 



that ca(*o) an d cb(<o) are not too close to 2 (recall that Ki are bounded away from 
zero and infinity, since (2.5) is n- variable mass-action system). The case when c^(io) 



is very close to 2 is not of interest to us as we want to show that ca cannot become 
too small, and therefore we look what happens when ca(£) is close to zero. 

The set C\ B is a positive translation of the 



Now we refer back to Figure 



4.2 



positive quadrant located at a small distance d from each of the axes, and with 
a cut at the corner near the origin. To illustrate the point, let us make a gross 
oversimplification and assume that C\ B is a square. (Note, however, that it is not, 
and, although the cut near the origin can be made arbitrarily small, it still requires a 
delicate analysis). With this simplification in place, we argue that ca cannot become 
smaller than d. Indeed, if at time t — t , the trajectory (ca, cb) reaches the boundary 
of C\ B and ca — d, then, as explained above, if c B (t ) is not too close to 2, the 
trajectory is pushed inside C\ B and ca increases. 

On the other hand, Theorem 4.1 does not apply for the projected system (5.13) 
if Cb(£o) is close to 2. However, in this case ccito) = 2 — c B (to) is small and we may 
project onto {A, C} instead. The projected reaction system 



O^A + C^A^C 2A^A^0, 
has rate constant functions 

K 1 (t) = K 1 (t)(2-C A (t)(2-C C (t)), K2{t) = K 2 {t), K 8 (i) = *«3(*), 
K 4 (t) =K 4 (t)(2-C C (t)), K 5 (t) = K 5 (t), K 6 (t) = K 6 (t)(2 ~ C A (t)) 2 

which arc all bounded away from zero at t = to. If C\ c C K^" 4 '^ is constructed in 
the same way as C\ B and at the same distance d from the coordinate axes, then, 
since c^(to) = d, we have (ca(*o)i c c(*o)) S ®£~a C anc ^ Theorem 
vector field at t — to points inside C\ c . Once again, ca must increase. 

One may recast the discussion above by using a symmetric construction of an 
invariant set T + . Namely, one considers the cylinder C\ 



implies that the 



A,B 



[c,d] anc j gjjniiar 

cylinders coming from all possible projections to pairs of variables. Defining T + to 
be their intersection, the reasoning above translates into T + being an invariant set 



for T(cq). This is precisely what we do in the proof of Theorem 5.1 



Note, however, that while the previous exposition sheds some light on the basic 



idea of the proof, (presented in section 5.4), the technical details involved are subtle 



and require an extensive preparation, which is the object of section 5.3 In particular, 
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the parameters required in the construction of the sets C + need to take into account 



the geometry of S(c$ ) and must be chosen carefully. Lemmas 5.1 and 5.2 are part of 
this process. 

5.3. Further preparation. As anticipated in the discussion above, a special 
case of Theorem |5.1| follows in a more or less straightforward way from Corollary |4.1| 

Proposition 5.1. Let (A/ - , k) = (S ,C,1Z, k) be a n-variable mass-action sys- 
tem with two-dimensional stoichiometric subspace and lower- endotactic stoichiomet- 
ric subnetworks. If the stoichiometric compatibility class S(cq) contains the origin, 
then T(cq) is a persistent trajectory. 

Proof. We denote W = {1, . . . , n}, so that the origin is the vertex Fw of S(co). 
Let (k,l) be a fixed pair in W 1 x W 2 , and let v*(W) = (pi, . . . ,p n ) and v 2 (W) = 
(qi , . . . , q n ). Note that all pi and may be chosen to be rational because S is generated 



by vectors of integer coordinates. As explained in Remark 5.2 we may assume that 
k < I and pi = qu = 1; note that we also have pt > 0, qi > for all i € {1, ... , n} 
from (5.1 1. 

Let A/" s = (S,C B ,1Z S ), s G {1, ... ,p} denote the stoichiometric subnetworks of A ". 
If \& denotes the matrix with columns v 1 (W) and v 2 (W), then \1/ has the form (4.1 ). 
As explained in Remark 5.3 for each s G {1, . . . ,p} we may choose a s e E™ such that 
C s cS + n s and nk,i(a s ) 



(0, 0). Since c(t) € S for all t > we have 



c(t) = *7Tfc,l(c(t)) 



(5.16) 



by (5.3) Remark 5.2 We have 



s=i p^p'en s 



and since, as implied by (5.6) Remark 5.4 the only reaction in 1Z S that projects onto 



Q — s- Q' via TTfc j is *Q + a s — >• ^PQ' + a s , it follows that 



3=1 Q^Q'e7r fc ,((K s ) 



where we denoted 



i s,Q- 



>Q'(t) = K*Q+a„ 



>*Q'+a. 



(*) 



for all Q — >• Q' e TTk,i(7Z s ). Therefore 7Tfc i ;(T(co)) is the trajectory of the 2D-reduced 
mass-action system [J P s=1 (^'k,i(J^'s), ^, k s , a s ) with initial condition irk,i{co). 

Corollary 



Since all 7r/ Ci ;(A r s ) are lower-endotactic by Proposition 3.1 



4.1 



implies 

that coordinates Xi(t) and Xk{t) of c(t) stay larger than some d > for all t > 0. 



Finally, (5.16) implies that all the coordinates of c(t) remain bounded away from 
zero. □ 

The special case of Theorem |5 . 1 1 contained in Proposition |5 . 1 1 illustrates well how 



projected systems, 2D-reduced mass-action systems and Theorem 4.1 come into play. 



As anticipated in section 5.2 the general case requires yet a little more preparation, 
which we discuss next. 

Fix a bounded trajectory T(cq) of (Af, k) and let M > 1 be such that T(c ) C 
[0, M] n . As hinted in section 5.2 we construct a polyhedron T + C M™ that stays 
away from d. 



™ and such that T(cq) C T + ; in the process we use the tools we have 
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developed thus far. Namely, we project (Af, n) onto well-chosen sets of two variables, 
we cast the projected system as a 2D-reduced mass-action system and we construct a 
corresponding C + set in each such two-dimensional face of R> - Finally, we construct 
certain cylinders out of the C + sets and define T + as the intersection of these cylinders. 

The projections to consider are of the form irk.i = it{k.i} with (k,l) G W 1 x W 2 
for all vertices Fw of S(cq)- 

We fix a vertex Fw and a pair (fc, I) G W 1 x W 2 . Based on Remark 5.2 if 
v 1 {W) — (pi, . . . ,p n ) and v 2 (W) — (qi,...,q n ) then we may assume that I < k 
and that pi = qk = 1- Note that if (k,l) G W 1 x W 2 then there is no other vertex 
F v of S(c a ) such that (k,l) G V 1 x V" 2 . Indeed, if F v = (/i, . . . , /„) then we have 
F v -F w = /jV 1 (W) + /fcV 2 (W). Since fc, I G Vi U F 2 C V, we have /; = f k = and 

SO Fy = . 

Recall that 7V S = (6>,C S ,7£ S ) denote the stoichiometric subnetworks of Af and 5 



denotes the stoichiometric subspace of Af. As in Remark 5.3 for each s G {1, .. . ,p}, 
let a s G E" such that C s Ca s + S. Let 

Q=(v 1 (W) v x (W)) 

be the matrix with columns v 1 (H / ) and v 2 (W) and define 

*fc,i = 7T W 6 = (tt w (v 2 (W)) tt w (v 2 (W))) (5.17) 

to be the matrix with columns 7Tvk(v 1 (W / )) and ir\y ^(v 2 ! (W)) . Note that and qi are 
non-negative for all i G and moreover, they are rational numbers since the stoi- 
chiometric subspace 5 of TV is generated by vectors of integer coordinates. Therefore 



^ k ,i is of the form (4.1) 



Since, by Proposition 3.1 7Tfc i(.A/" s ) is lower-endotactic for all s G {1, . . . ,p}, we 
may construct the set 

= ^(W(^ s )}i<,<„,*M>V,M) c Z C{M} (5.18) 



as in section 



4.3 



such that nf. z(co) G ^jtr ^e will choose ?y' in what follows; also, we 
will take advantage of the flexibility in the construction of C\ l to equip this set with 
a few useful technical properties. 

We start with two lemmas which show the intuitively clear facts that if a point 
of S(cq) is close to Zyy then it is also close to Fw, and that if some components of a 
point in S(cq) are small, then the point is close to a face where all those components 
are zero. 

Lemma 5.1. Let Fw be a face of S(c$). Then there exists 6 > such that 
dist(x, Fw) < #dist(x, Zw) for all x G S{cq). 

Proof. If S(cq) has dimension one, we denote by a the angle between S(cq) 
and Zw ■ Since S(cq) intersects the positive orthant we have a > and therefore 
< sec a < oo. Then dist(x, Fw) — <5dist(x, Zw) for all x G S(cq), where S = sec a. If 
S(cq) has dimension two, for x G S(cq) we let pf and pz denote the projections of x on 
Fw and Zw- Let a = inf xg s( Co ) Z(xp^p^). If a = 0, then, since S(cq) is closed, there 
exists xo G S(cq)\Fw such that Z(xoP_fPz) = 0. In turn, this implies that pp = pz 
and therefore S(cq) contains the line segment XoP_f which is perpendicular to Zw- It 
follows that dist(x, Fw) = dist(x, Zw) for all x G S(cq). If a > then we let 5 = sec a. 
For x G 5(co) we have dist(x, Fw) = sec Z(xpi?pz)dist(x, Zw) < <5dist(x, Zw) for 
all x G S(c ). □ 

Lemma 5.2. There exists A > such that if I C {1, . . . , n} and d = (d\, . . . , d n ) G 
S(co) is such that di < X for i G I, then for some face Fw of S(cq) we have I C W. 
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Proof. If the origin is a face of S(cq) the claim in the lemma is clearly true 
(any positive value for A will do). Otherwise, for J C {1, . . . , n} we define m(J) = 
inf xS 5( CQ ) dist(x, Zj) and we let A = min{m(J) | J C {1, . . . ,n},m(J) > 0}/y/n. We 
have A > and 



m(I) 2 < dist(d, Zjf = ^ dj < n\ 2 



iei 

which shows that m(J) = and the conclusion follows. □ 

In view of Proposition 5.1 we may assume that the origin is not a vertex of S(c ). 
Let v mm denote the smallest nonzero coordinate of a vertex of S(co) and, fixing a A 
given by Lemma [B~2"l let 

C = min{A,t> mm /2,l}. (5.19) 
Moreover, let 1 = (1, .... 1) G E", E = max (P ■ 1) and define 



7] = mm 



in{^,^}. (5.20) 



Recall from section 4.3 that the construction of a set C + depends on the numbers £ 
and d. Also recall that £ may be chosen arbitrarily small; once £ is fixed, d may also 
be made small independently of the value of £. Since there are finitely many pairs 

{k,l) G W 1 x W 2 (counting all vertices Fy/ of S(cq)) we can choose the same values 

+ 

A:,r 



of £ and d in the construction of all sets We fix £ small enough such that 



|J {\ Pi / Pj \ I z, j G Cl^ 1 } U | i, j G CW 2 } k < 

{.Fw vertex of S(cq) ) 
v 1 (W)=(pi,...,p„) 
v 2 (W) = ( 9 i,...,g„)} 



4 



(5.21) 

As can be seen from Figure 4.2 the shape of £k,i near (0, 0) enables us to choose 
e > such that 

[0, e] 2 n £+, = for all vertices F w of S(c ) and all (fc, I) e W 1 x W 2 . (5.22) 
We now choose d such that 

d<min|A,C,e (J ({ft/Pj I <»J e W\W 1 } U | i,j G W\W 2 }) I, 

l {i*V vertex of S(co) / 

v 1 (W) = (p 1 ,...,p„) 
v 2 (W) = ( 9 i,...,g„)} 

(5.23) 

where we recall that A is given by Lemma 5.2 £ is defined in (5.19) and e was chosen 
to satisfy ( |5.22[ ). 

We shift C K i (and £+,) along R c i k ^ and define 

H fe>i = {x G R n | (sfc.a:,) G £ M } and «+ = {x G R" | (a*,^) G £+ }. 

Note that "H^ z = conv('Hfc.z)- Finally, let 

T+ = (rfl+M| )n P| f| U+ v (5.24) 

F w vertex of S(c ) (k,l)£ W 1 X W 2 
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By definition, the convex polyhedron T + lies in a positive translation of the nonneg- 
ative orthant. In view of Theorem |2.1[ we shall be concerned with the behavior of 
the flow c(t) on the boundary of T + ; we conclude the preparatory discussion with 
the following lemma, which shows that the part of <9T + that is of interest to us 
does not include the boundary of (dl + R> )j but only the boundaries of H~£ z , for 

{k, I) e w 1 x w 2 . 
Lemma 5.3. 

dT + nS(c Q )c |J |J (%,,nS(c )). 

F w vertex of S(c ) (k,l)eW 1 xW 2 



Proof. We have 

n 

dT+C \J(dl + Z {t} )U 



U U KM- 

F w vertex of S(c ) (k,l)eW 1 xW 2 



Suppose x e <9T n S'(co) and x £ dl + ^{n- Since d < X (see (5.23)), Lemma 5.2 



implies that there exists a face F\y such that i £ W. Possibly making W larger, we 
can assume that Fw is a vertex of S(cq). Now we show that i £ W 1 U W 2 ; suppose 
this was false and let (k,l) £ W 1 x W 2 . Assume that I < k (otherwise swap W 1 with 
W 2 ) and let w 1 {W) — (pi, . . . ,p n ) and v 2 (W) = (qi, . . . , q n ) where pi = qk = 1- S ince 
i £ W^W 1 U W 2 ), both pi and qi are strictly positive, as explained in Remark 
Using (5.23) we obtain 



5.1 



uiax{x ipi,x k q t } < xipi + x k qi = Xi = d < min{(q i /q fc )e, (j)i/pi)e} = mi^^e,^^} 



and so (xi,Xk) £ [0,e] . This, together with (5.22) implies that x ^ an< ^ therefore 



x ^ T + , contradiction. We conclude that i £ W 1 U W 2 . Suppose i £ W 1 , and let 
j £ W 2 . We have 

x6K+.n(c(l + %) CH id . 



5.4. Putting things together. We are ready to prove our main persistence 
result, Theorem |5.1| We keep using the notations introduced thus far in this section; 
in particular, if Fw is a vertex of S(co) and (fc, I) £ W 1 x W 2 we assume that I < k, 
we denote v 1 (W) = (pi, . . . ,p n ),v 2 (W) = (qi,. . . ,q n ) and we set pi = q k = 1. 



Proof of Theorem 5.1 For any vertex Fw of S(cq) and for any pair (k,l) £ 
W 1 x W 2 we have by construction c G Hti> an< ^ therefore c € T + ■ We will use 
Theorem 2.1 to show that T(cq) C T + ■ Suppose t > is such that x = c(t) £ dT + ; 



it is enough to show that 



n ■ c(t) > 



for all n £ —N-j-+ (x). Lemma 5.3 implies the existence of a vertex Fw = 



k,l- 



(5.25) 

■ ■ i fn) 

The 



of S(cq) and the existence of a pair (fc, Z) € W 1 x TU 2 such that x £ H 
generators of the convex cone — A^7-+(x) lie in the union of —Nn,+ (x) for all pairs 



(k,l) such that x £ %k,i, thus (5.25) needs only be verified for vectors n belonging to 
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this union. Therefore we fix such a pair (k,l) and we let n 6 —N^+ (x). Since only 
the coordinates k and I of n are nonzero, inequality (5.25) is equivalent to 



7Tfc,*(n) • 7Tfc,;(c(t)) > 0. 



(5.26) 



According to Remark 5.3 for each s e {1, .. . ,p} we may choose a s € R ra such that 
C s cS + n s and nk,i(a s ) = (0, 0). 

Case (i). Suppose 7Tfc /(x) lies on the finite part of For c € T(cq) we have 
c — € S and we may write, using Remark 5.2 and the definition (5.17) of ^k,i ■ 

n w (c) = tt w (c - F w ) = ir w (QTr k j(c- F w )) = $ ktl n k j(c - F w ) = * feii 7r fei /(c). 

Therefore ttw{c) — ^k,i^k,i{c) for any c € T(co). On the other hand, according 
to Remark 



5.4 



the only reaction in TZ S mapped by nk,i to Q — > Q' € itk,i(JZ-s) is 
QQ + a s — > QQ' + a s . It follows that the dynamics of (AT,k) projected onto {k,l} 
may be written 



dt 
p 



tm(c)=E E *P-^P'(*)'n!w(c) ,rBwr(i 'W(c) ,rv(P) (7r fc , I (i y )-7r fc , I (P)) 

s=l P^P>£K S 



E 

8=1 I 

E 



£ K e Q+ a^eQ'+a s (t)7r Cw (c)^-( 0c 3+^)7r w (c)^'( e ^)(Q' - Q) 
2'eTTfc,i(7i,) 



8=1 I 

where 



£ K s , Q ^ Q ,(<)(vI/ fc , 7 r fc , i ( C ))' J/fe ' !C3+a =(Q' - Q) 



(5.27) 



*«,«-»<J'(*) = «eQ + a^e Q ' + a s (i)7r CM ,(c(t))^«- (0Q+a ^ for all Q -> Q' e 7T fe ,/(^ s ). 

(5.28) 

Therefore the system of differential equations (5.27) is the 2D-reduced mass-action 
system \J P s=1 (nk,i(Af s ), #k,i, ««,a a ). 

For all i G CW 7 we have (recall (5.4)): 

Xi = ft + xipi + x k qi > fi- \pi\xi - \qi\x k = fi- \pi/pi\xi - \qi/qk\xk- 

Moreover, since 7Tj, ;(x) lies on the finite part of Ck,\ we have Xk, xi < £. This, together 
with (5.21 ) implies, for any i € CW : 

ft - \Pi/pi\xi - \qt/qk\x k > fi - (\Pi/pi\ + Ift/fffcDe > fi ~ v mm /2 > V min /2. 

Therefore, for all i € ZW we have 

Xi>fi- v mm /2 > v mm /2 

(recal l that v mm denotes the minimum of nonzero coordinates of Fw)- It then follows 
from (5.19) that X{ > £ for any i 6 CW. This yields 



M (BQ+a s 



> 



( x )^Cw(eQ+a s ) > C (6Q+a s ).l for &ny Q ^ Q' £ 7^,(^)5 



recalling (5.28) and using (5.20), we then have k s ,q->-Q' (i) G (Vi V 7 /)- Sin ce 7T fc,;(n) € 



-#r+/^MW) and £+, = £+ ; ({7r fcii (A4)}i<.< P , * fe ,i. V, M )> Theorem 
( pT26| '.' 
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4.1 



implies 



Case (ii). Now suppose that 7rfc,j(x) lies on the infinite part of C^^y, for instance 
Xf. = d. Let I = {i € {1, . . . ,n} \ Xi < Q an d note that, by (5.23) we have k E I. 
Since, by (5.19) we have £ < A, Lemma 5.3 implies that there exists a face Fgy of 

! | 2 

S(cq ) with / C W . We may assume that F^r is a vertex. We claim that k & W L)W ; 

indeed, otherwise, let (k, I) G W x W and v 1 (W) = (pi . . .p„), \ 2 (W) = (qx . . . q n ) 
such that pi = gj = 1. Since, by (5.22), at least one of and is larger than e, in 



view of (5.23) we have the following contradiction: 



d = x k = xfPh +xj,q k > emin{p k ,q k } > d. 



Therefore k S W U W ; suppose k € W and let 7 G W . Since for each i e CW 
we have > C (this from our definition of /), the same argument as in case (i) shows 
that 7r fe j(n) • 7r fc r(c(t)) > 0. On the other hand, since only the fc-th coordinate of n is 
nonzero, we have 



7rfc,;( n ) ' T k ,l(.c(t)) = vr, |(n) • n k j(c(t)) 



and (5.26) is shown. 



Recall that weakly reversible reaction networks are endotactic and in particular 
lower-endotactic. The following corollary states that the version of the Persistence 
Conjecture proposed in [2] holds for systems with two-dimensional stoichiometric sub- 
space. 

Corollary 5.1. Any bounded trajectory of a weakly reversible n-variable mass- 
action system with two-dimensional stoichiometric subspace is persistent. 

Example 8. To conclude this section let us revisit the n-variable mass-action 
system (2.5). We know that its stoichiometric subspace is two-dimensional (Example 



1) and that its stoichiometric subnetworks coincide with its linkage classes (Example 
4)- If L± denotes the first linkage class then the projection 71^2 : aff(Li) — > -^{3,4} 
is invertible. Since 7ri a( Li) = {B ^ A — > A + B — > 0} is easily seen to be endotactic, 
according to Proposition 3.1, the same is true for L 1 = 2 { 7T l,2(Ll)) ■ Similarly, 
the second linkage class is endotactic. Since ca + cd and cb + cq remain constant 
along trajectories, any trajectory is bounded. Therefore, Theorem\5.1\ implies that the 



dynamical system (2.6) is persistent 



6. The Global Attractor Conjecture for systems with three-dimensional 
stoichiometric subspace. Recall from Introduction that in order to show the Global 
Attractor Conjecture it is enough to prove that all trajectories of complex-balanced 
mass-action systems are persistent. Theorem |5.1| may be used to analyze the behavior 
of trajectories of weakly reversible mass-action systems near faces of S(cq) of codi- 
mension two. As we shall see below, trajectories can approach such a face only if they 
approach its boundary. This is made precise in Theorem |6.2| On the other hand, as 
discussed in Introduction, vertices of S(cq) cannot be cj-limit points for trajectories 
of complex- balanced systems [TJ [5], Moreover, codimension-one faces of S(cq) are 
repelling [1] and we have the following result: 

Theorem 6.1 (@] Corollary 3.3). Let c e M" and let T(c Q ) denote a bounded 
trajectory of a weakly reversible complex-balanced mass-action system. Also let Fyy be 
a codimension-one face of S(cq). IfT(cg) does not have u-limit points on the (relative) 
boundary of F\y, then it does not have cu-limit points on Fyy. 

For weakly reversible, complex-balanced systems with three-dimensional stoichio- 
metric compatibility classes, the results mentioned above cover all faces of S(cq) and 
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can be combined into a proof of the Global Attractor Conjecture for this case. We 
start with the following lemma. 

Lemma 6.1. Let (S,C,1Z, re) be a weakly-reversible n-variable mass-action system, 
let Co € M>o Q> n d let F\V oe a f ace °f S(cq) of codimension two. Then for any compact 
K C int(_FV) there exist r > and e > such that if for some t! ,t" £ R>o we have 
c(i) = ( Xl (t),...,x n (t)) € tt Cw (K) x [0,e] w for all t G [t',t"], then E ieH /^( i ") > 
T ^2iew x iV" )• 

Proof. We denote Af = (S,C,1Z). Let S denote the stoichiometric subspace of 
Af, let d = dim 5*, and let irw\s '■ S — > Z$ w denote the restriction of ttw to S. 
Since Fw is of codimension two, we have dim(5 n Zw) = d — 2. But S n Z\y = 
ker(ir\y\s) and therefore ttw(S) = Im(irw\s) has dimension two. Note that ttw(S) 
is the stoichiometric subspace of A^ = ww(Af). Since Af is weakly reversible, so is Af; 
in particular, Af has two-dimensional stoichiometric subspace and lower- endot act ic 
subnetworks, which we denote by Af s = (S,C a ,1Z s ) for s € {1, . . . ,p}. 

We may now apply the results of the preceding sections. The face Fyy of the 
stoichiometric compatibility class S(nw(co)) of Af is the origin of R w . We let F w i 
and F W 2 denote the two edges of S(nw(co)) that are adjacent to F\y (recall that 
W 1 , W 2 are subsets of {1, ... , n} which are contained in W). Let (k,l) € W 1 x W 2 . 
We may assume I < k and we scale the direction vectors vi(W) and v 2 (T / K) of F w i 
and Fy^i such that the Zth coordinate of Vi(W) and the fcth coordinate of V2(W) are 
equal to 1. Also, for each s e {1, . . . ,p}, we choose a s £ R w such that C s C a 8 + S 
and TTk,i(a s ) — (0,0) (this is possible as explained in Remark 5.3). 

Suppose that re(t) £ (rj, l/i]) 71 for all t > and let ^k,i be the matrix with 
columns v\(W) and ^(W). We have 



|tt m (c)= E K P ^ P .{t)c p {Q' -Q) 

7Tk,l(P-+P') = Q^Q'} 

E E E n P ^ P ,(t)c p (Q'-Q), 

s=l Q^Q'e-n k .i(n 3 ) {P^P'£1Z\ 

Trw{P^P')=Vk,iQ+a 3 ^9k,iQ'+a a } 



since, as explained in Remark 5.4 the only reaction of 1Z S that is mapped to Q — > Q' 
by 7Tfe,i is $>k,lQ + a s *' fe; ;Q + a s . Further, note that tt w (c) = ^k,l^k,l( c ) for am/ 
c e T(co); therefore, writing c p = 7r 1 v(c) 7rw( - p '7r[; vl/ (c) 7rc «'(- p ) we have 

j t ^Ac) = jz E ^, Q ^W(*M^(c))* fc ' !Q+as (Q'-Q), 



where 



H 



s,Q- 



'(*)= E K p^(i)7T CH/C (t)^- 



{p^p'eK| 

7r w (P^P')=*fe,iQ+«a^*fe, ( Q'+«»} 



for all t > 0. Since if C int(Fw) there exists < £ < 1 such that ir$w(K) C 
(C, l/C) Cw . It follows that there exists 77' < 1 such that re S) g_>.Q/(t) € (?/, l/V) for all 
£ G The projection of (Af, re) onto {fc,Z} is a 2 D-red uced mass-action system 

Ug =1 (7Tfe > i(A/' s ),Vl/fc ! ;,re s ,a s ). Since fc, Z € W 7 , Corollary 4.2 implies that there exists 
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e > and t' > such that if c(t) E tt Cw (K) x [0, e] w for t € [t',t"}, then x k (t") + 
xi{t") > T'(x k (t') +X[(t')). Since n w (c(t)) = * fe ,;7r fci/ (c(i)) we have £ ieVF £i(i) = 
Aajfe(t) + Bxi(t) for some positive numbers A and £? and for all i > 0. Therefore 

V Xi (t") > mm{A,B}(x k (t")+ Xl (t")) > r'min{A,B}(x k (t')+xi(t')) > T ' mm { A ^> V xtf) 
— ' max A, B \ c — ' 

and the conclusion follows by setting r = ~~rj^gr^ ■ 



We may now extend the result of Corollary 3.3 in [4] (Theorem 6.1) to faces of 
codimension two. 

Theorem 6.2. Let (S,C,1Z, k) be a weakly-reversible n-variable mass-action sys- 
tem, let Co € M" suc h that the forward trajectory T(cq) is bounded, and let Fw be a 
face of S(cq) of codimension two. Then, ifT(co) has uj -limit points on Fw, it must 
also have u -limit points on the relative boundary of Fw . 

Proof. Suppose the claim of the theorem was false. Then lim w T{c$) n Fw C 
int(i*V) is compact (it is an intersection of closed sets and is bounded since T(cq) 
is bounded). It follows that there exists a compact set K C int(i'V) such that 
lim w T(cq) n Fw C int(if). From Lemma [6. 1| there exist r > and e > such that 



if c(t) = (xi(t),...,a? n (t)) € vr CH/ (X) x [O^T for all * e [f,f] then E^a^t") > 
T J2iew x i(t')- It follows that, in order to approach an w- limit point in K, T(cq) 
must exit and reenter K e = n^(K) x [0,e] infinitely often. More precisely, there 
exist t\ < t% < ... with t m — > oo as m — > oo such that c(t rn ) G dK t for all 
m > 1 and E»ew a '»(*m) as m -) oo. Note that K = ^^(Jsr) x {0} w 
and therefore we must have c(t m ) £ d(n^ w (K)) x [0, e] w for m large enough. If 
A = {c(ti), c(<2), . . .} then, since the closure A" of A' in K ra is compact, there must 
exist y € A C ^(^^(if)) x [0,6]^ such that Eiew^ = 0, which implies y € 
5('rcw(- fs: )) x {O} 1 *' = <9if (the relative boundary of if). But A n F w = since 
T(co) does not intersect the boundary of M. n , and it follows that y is an accumu- 
lation point of A, and therefore y G lim w T(co). But this is a contradiction since 
lirriu T(c ) n d# = 0. □ 

A proof of the Global Attractor Conjecture may now be obtained for systems 
with three-dimensional stoichiometric subspace. 

Theorem 6.3. Consider a complex-balanced system with three-dimensional stoi- 
chiometric subspace. Then, the unique positive equilibrium contained in a stoichiomet- 
ric compatibility class is a global attractor of the relative interior of that stoichiometric 
compatibility class. 

Proof. It is known that any trajectory of a complex-balanced system is bounded 
[T2] ■ Since a stoichiometric compatibility class has only faces of codimension two, 
codimension one and vertices, Theorems |6 . 1 1 and |6 . 2[ applied in this order, show that 
if a trajectory T(co) has w-limit points, then a vertex of S(cq) is an w-limit point. 
But this is known to be false [U [9] . □ 
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